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Section  1 
INTRODUCTION 

The  basic  problem  addressed  in  this  report  is  the  uncertainty  in 
measured  stagnation  enthalpy  in  the  Reentry  Nose  Tip  (RENT)  Leg  of  the 
AFFDL  50  Megawatt  Facilities.  The  values  of  local  stagnation  enthalpy 
derived  from  stagnation  point  heat  transfer  measurements  appear  to  be 
significantly  higher  than  those  obtained  from  heat  balance  measurements 
on  the  arc  heater  and  from  direct  measurements  of  local  stagnation  en- 
thalpy. These  discrepancies  may  be  due  to  possible  errors  in  instrumenta- 
tion or  due  to  some  of  the  assumptions  contained  in  stagnation-point 
heat  transfer  theories  being  violated.  The  basic  objective  of  this 
program  was  to  investigate  the  latter  possibility,  namely,  the  possibility 
that  there  are  some  important  physical  phenomena  present  in  the  RENT,  such 
as  enthalpy  gradients,  free-stream  turbulence,  etc,,  which  are  not  taken 
into  account  by  the  Fay  and  Riddell  type  of  stagnation-point  heat  transfer 
theories. 

The  technical  approach  to  this  problem  was  to  critically  review 
stagnation  point  heat  transfer  theories  and  other  papers  and  reports 
which  may  be  relevant  to  a possible  cause  of  measured  discrepancies  in 
stagnation  enthalpy.  Two  groups  of  papers  and  reports  were  considered. 

The  first  group  consisted  of  stagnation  point  heat  transfer  theories  defined 
in  ?.  narrow  sense,  i.e.,  theories  treating  the  heat  transfer  from  a dissocia- 
ted gas  to  blunt  bodies  in  a supersonic  stream  (e.g.,  Fay  and  Riddell).  The 
second  group  consisted  of  papers  and  reports  which  deal  with  phenomena  such 
as  particle  impingement,  free  stream  turbulence,  etc.,  which  are  excluded  from 
stagnation  point  heat  transfer  theories  as  defined  in  the  above  narrow  sense, 
but  which  may  be  present  in  the  RENT  facility  and  may  be  a possible  cause 
of  the  observed  discrepancies  in  local  stagnation  enthalpy  measurements. 

This  second  group  is  referred  to  as  "possible  cause"  theories.  Also 
included  in  the  "possible  cause"  group  are  reports  and  papers  containing 
data  and  information  on  facility  operating  conditions  and  facility  test 
data  which  are  useful  in  evaluating  the  relevancy  of  possible  cause 
phenomena.  The  technical  approach  included,  in  addition  to  that  described 
above,  analyses  to  determine  the  relevancy  of  a possible  cause  and  investiga- 
tions to  determine  cause  and  effect  relationship.  Analytical  investigations 
of  the  possible  effects  of  timewise  fluctuation  of  flow  properties,  radial 


-1- 


enthalpy  gradients  and  swirl  flow  on  heat  transfer  were  conducted.  A cause 
and  effect  relationship  between  free  stream  turbulence  and  heat  transfer 


enhancement  v;as  established. 

A total  of  eight  possible  causes  were  investigated  in  depth.  The 
findings  from  these  investigations  are  summarized  in  Table  1-1.  (Radiative 
heat  transfer  was  also  briefly  investigated  and  dismissed  as  being 
negligible).  The  major  conclusions  of  these  investigations  are  that,  1) 
the  discrepancy  between  the  arc  chamber  heat  balance  enthalpy  and  derived 
free  stream  enthalpy  is  between  20%  and  50%,  and  2)  the  possible  cause 
which  may  be  responsible  is  free  stream  turbulence. 

Recommendations  for  additional  efforts  to  further  investigate  this 
problem  are  two-fold.  The  first  recommendation  is  to  experimentally  deter- 
mine the  free  stream  turbulence  level  and  the  copper  mass  loading  in  the  RENT 
facility.  These  flow  properties  are  critical  to  the  quantitative  evaluation 
of  potentially  important  possible  causes.  The  second  recommendation  is  to 
develop  two  theoretical  models.  The  first  model  would  predict  the  flow 


properties  in  the  test  section  (free  stream  turbulence,  total  pressure, 
enthalpy,  etc.),  including  the  effects  of  swirl  in  the  arc  chamber.  The 
second  model  is  a turbulent  viscous,  shock  layer  model  which  will  predict 
the  heat  transfer  enhancement  and  change  in  heat  transfer  distribution  due 


to  free  stream  turbulence  by  properly  accounting  for  the  interaction  between 


free  stream  turbulence  and  the  turbulence  generated  in  the  boundary  layer. 
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Section  7. 

STAGNATION  POINT  HEAT  TRANSFER  THEORIES 


In  this  section,  conventional  stagnation  point  heat  transfer  theories 
and  heat  transfer  theories  which  attempt  to  account  for  the  effect  of  free- 
stream  turbulence  are  described.  Although  the  examination  of  conventional 
heat  transfer  theories  has  contributed  little  to  the  understanding  of  the 
RENT  stream  enthalpy  discrepancy,  a brief  summary  of  these  theories  are 
included  for  the  sake  of  completeness. 


2 . 1 Conventional  Stagnation  Point  Heat  Transfer  Theories 


Conventional  stagnation  point  heat  transfer  theories  are  defined  as 
those  theories  which  predict  the  heat  transfer  at  the  stagnation  point  of 
a body  immersed  in  a uniform  flow  which  is  free  of  gradients  of  any  type, 
has  no  free  stream  turbulence,  and  is  devoid  of  any  type  of  particulate 
matter.  These  theories  have  been  developed  over  the  past  20  years  and 
are  generally  valid  for  supersonic  and  hypersonic  flight  velocities, 
dissociated  and  ionized  gases  and  for  flows  that  are  in  or  out  of  chemical 
equilibrium. 


In  this  section,  only  a general  description  of  stagnation  point  heat 
transfer  theories  are  presented.  A detailed  discussion  of  chemical  non- 
equilibrium effects  is  presented  in  Section  8 and  the  effects  of  transport 
properties  variations  on  heat  transfer  are  discussed  in  Section  9. 

Conventional  stagnation  point  heat  transfer  theories  which  are  of 
interest  to  the  present  program  are  typified  by  Fay  and  Riddell  (Ref.  2-1), 


Ref.  (2-1)  Fay,  J.  A.,  and  Riddell,  F.  R.,  "Theory  of  Stagnation  Point 

Heat  Transfer  in  Dissociated  Air",  J.  Aeronaut.  Sci.,  Vol.  25, 
No.  2,  Feb,  1985,  pp.  7T-85,  121 


I 

f 

I 
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Hoshizaki  (Ref.  2-2)  and  others  (Refs.  2-3  to  2-5).  In  all  of  these 


theoretical  treatments  of  stagnation  point  heat  transfer,  it  is  implicitly 
assumed  that  the  free  stream  is  uniform  and  turbulence  is  negligible.  The 
oncoming  flow  is  assumed  to  be  supersonic  or  hypersonic  so  tha  . a bow 
shock  is  formed.  The  Reynolds  number  is  assumed  to  be  sufficiently  high 
' that  the  stagnation  point  boundary  layer  thickness  is  much  less  than  the 

; shock  detachment  distance.  Under  these  conditions,  the  vorticity  generated 

j by  the  curved  shock  can  be  neglected  and  the  flow  at  the  edge  of  the 

I 

I boundary  layer  can  be  assumed  to  be  free  of  gradients  normal  to  the  surface. 

t 

I 

The  high  flight  velocities  considered  generate  stagnation  temperatures 


large  enough  to  dissociate  and  ionize  the  gas.  The  gas  is  assumed  to  be  a 
mixture  of  perfect  gases  and  in  most  treatments  it  is  assumed  to  be  either 
a binary  mixture  or  a multicomponent  mixture  with  equal  mass  diffusion 
properties.  Transport  properties  such  as  viscosity  and  thermal  conductivity 


are  assumed  to  vary  across  the  boundary  layer. 

At  the  stagnation  point,  the  flow  is  locally  similar  and  partial  dif- 
ferential equations  for  the  conservation  of  mass  (overall  and  species) 
momentum  and  energy  can  be  reduced  to  a set  of  total  differential  equations 
with  the  similarity  variable  as  the  independent  variable.  Numerical  solu- 
tions to  tliese  equations  require  iteration  on  either  the  initial  conditions 


Ref.  (2-2)  Hoshizaki,  H.,  ''Heat  Transfer  in  Planetary  Atmospheres  at  Super- 
Satellite  Speeds",  ARS  J.,  Vol.  32,  No.  10,  Oct.  1962,  pp.  1544- 
1552 

Ref.  (2-3)  Pallone,  Adrian  and  Van  Tassell,  Wm. , "Effects  of  Ionization  on 
Stagnation  Point  Heat  Tra >sfer  in  Air  and  in  Nitrogen",  Phys. 
Fluids,  Vol.  6,  No.  7,  July  1963,  pp.  983-986 
Ref.  (2-4)  Marvin,  J.  G.,  and  Deiwert,  G.  S.,  "Convective  Heat  Transfer  in 
Planetary  Gases",  NASA  TR  R-224,  1965 


at  the  wall  (for  integrating  from  the  wall  to  infinity,  in  similarity 
coordinate)  or  on  entire  profiles  if  the  equations  are  cast  in  integral 
form  (method  of  successive  approximations). 


The  transport  of  chemical  enthalpy  can  be  handled  in  two  ways  in  these 
solutions.  One  can  retain  the  terms  which  account  for  chemical  energy 
transport  in  the  energy  equation  and  solve  the  equations  in  this  form  (Fay 
and  Riddell,  Ref.  2-1),  or  one  can  define  a total  thermal  conductivity 
which  is  the  sum  of  the  thermal  conductivity  due  to  molecular  collisions 
and  the  thermal  conductivity  due  to  diffusion  of  atoms  or  ions  (Koshizaki, 

Ref.  2-2).  The  latter  method  is  restricted  to  thermochemicai  equilibrium 
and  a fixed  elemental  gas  composition. 

It  has  been  a common  practice  to  correlate  the  numerical  results  v?ith 
total  enthalpy  or  flight  velocity,  fluid  properties  evaluated  at  the  wall, 
or  stagnation  pressure.  The  correlations  that  have  been  published  are  all 
quite  similar  in  form  and  all  yield  stagnation  point  heat  transfer  predic- 
tions which  vary  by  about  5 to  10  percent.  These  differences  are  usually 
attributed  to  slight  differences  in  transport  properties  and  numerical 
methods  of  solution. 

Analysis  of  these  theories  leads  to  the  conclusion  that  the  RENT  enthalpy 
discrepancy  cannot  be  explained  on  the  basis  of  gross  inacoiuctLius  in  the 
theories.  The  theories  have  been  verified  by  comparison  to  iabciatory  and 
flight  test  data  and  have  been  employed  in  the  analysis  and  development  of 
many  successful  reentry  vehicles.  The  theories  are  valid  within  the  frame- 
work of  the  assumptions  and  restrictions  for  which  they  were  developed, 

2,2  Stagnation  Point  Heat  Transfer  Theories  with  Free  Stream  Turbulence 

The  stagnation  point  heat  transfer  enhancement  due  to  free  stream  tur- 
bulence has  been  studied  theoretically  for  incompressible  flow  by  Galloway 
(Ref,  2-6)  and  Traci  and  Wilcox  (Ref,  2-7).  Galloway  includes  a Reynold.^ 


Ref.  (2-6)  T.  R.  Galloway,  "Enhancement  of  Stagnation  Flow  Heat  and  Mass 

Transfer  Through  Interactions  of  Free  Stream  Turbulence",  AlChE 
Journal,  Vol,  19,  No.  3,  1973 

Ref,  (2-7)  R.  M.  Traci  and  D,  C.  Wilcox,  "Analytical  Study  of  Freestream  Tur- 
bulence Effects  on  Stagnation  Point  Flow  and  Heat  Transfer",  AIAA 

7th  Fluid  and  Plasma  Dynamics  Conf.,  Palo  Alto,  Ca.,  June  17-19, 
1974,  AIAA  Paper  No.  74-515 
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stress  term  in  the  otherwise  I'uinar  boundary  layer  questions.  This  Rey- 
nolds stress  was  modeled  by  a somewhat  modified  mixing-length  eddy  viscosity 
turbulence  model.  It  was  assumed  by  Galloway  that  the  velocity  fluctuations 


!vl  = 


by 


(2-1) 


where  i is  the  Prandtl's  mixing  length  = ky  and  k is  a constant. 


It  was  also  assumed  that  the  Reynolds  stress  u'v'  is  proportional  to 


V 


= u- 


, bu  bu 


(2-2) 


where  c is  defined  as  the  eddy  viscosity.  'Jnlike  the  conventional 
m 

mixing  length  turbulence  model,  which,  according  to  Eqs.  (2-1)  and 

(2-2),  should  be  given  as 

_ , 2 2 bu. 
e = k y (— ), 
m by 


(2-3) 


Galloway  created  a relationship  between  and  the  free  stream 

turbulence  as 


e = k lU  y (2-4) 

where  I is  the  free  stream  turbulence  intensity  and  U is  the  free 

CO 

stream  velocity.  It  is  to  be  noted  that  should  Eq.  (2-3)  be  used,  there 
would  be  no  relation  between  boundary  layer  Reynolds  stress  and  the  free 
stream  turbulence.  Also  there  would  be  no  stagnation  heat  transfer 
enhancement  due  to  the  fact  that  u = 0,  = 0 at  the  stagnation  region. 

Hence,  the  conventional  mixing  length  model  is  not  an  applicable  tur- 
bulence model  for  the  study  of  free  stream  turbulence. 

Based  on  the  eddy  viscosity  of  Eq.  (2-4),  enhancement  in  heat  transfer 

at  the  stagnation  point  was  found  by  Galloway  as  shown  in  Figure  1.  In  this 

figure,  the  temperature  profiles  normal  to  the  surface  are  shown  as  a 

function  of  the  parameter  X which  is  equal  to  kx7R  , where  R is  the 

e e 

Reynolds  number  based  on  nose  radius  R.  The  Nusselt  number  is  also 
tabulated  in  Figure  1, 
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Tj  = Y/R  • 


Figure  1 Profiles  Normal  to  the  Surface 
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Comparisons  with  experiniental  results  were  also  made  by  Galloway.  By 
choosing  the  right  value  of  the  constant  k (Eq . 2-2),  good  agreement  between 
theory  and  experiments  were  found  as  discussed  in  Section  4 (Free  Stream 
Turbulence)  of  this  report. 

Based  on  the  work  of  Ref.  (2-6),  we  may  conclude  that  free  stream  tur- 
bulence significantly  enhances  the  heat  transfer,  the  relevant  parameter 
for  this  enhancement  is 

Traci  and  Wilcox  (Ref,  2-7),  on  the  other  hand,  used  a turbulent 
kinetic  energy  (TKE)  model  in  their  study.  The  TKE  model  they  used  was 
developed  by  Saffman  (Ref,  2-8).  In  this  model  the  relation  between  Reynolds 
stress  and  the  mean  rate  of  strain  tensor  is  given  as 


7..  - 2(v  + e ^ S..  - ■re  0., 
ij  m ij  3 ij 


(2-5) 


where  t..  is  the  Reynolds  stress  tensors.  S..  is  the  mean  rate  of  strain 
ij  ij 

tensor  which  is  related  to  the  gradients  of  the  mean  flow,  is  the 
Kronecker  delta  (i,e.,  b^^  = 0,  if  i?^j,  b^^  = 1 if  i = j),  v is  the 
molecular  viscosity,  is  the  eddy  viscosity  and  e is  the  turbulent 
kinetic  energy  density. 

e is  related  to  the  turbulent  kinetic  energy  e and  the  so-called 
m 

turbulent  pseudovorticity,  tr,  as 


e 

^m  tu 


The  turbulent  pseudovorticity  can  be  related  to  the  turbulent  length 
scale  as 


(2-6) 


(2-7) 


where  a*  's  a constant. 

To  complete  the  model,  partial  differential  equations  for  e and 
it  are  derived.  These  equations  state  the  conservation  of  turbulent  kinetic 
energy  and  turbulent  pseudovorticity  along  a streamline.  The  solution  is 
to  be  sought  by  coupling  these  equations  with  the  flow  equations  , 


Ref.  (2-8)  P.  G.  Saffman,  "A  Model  for  Inhomogeneous  Turbulent  Flows, 
Proceedings  of  Royal  Society,  A317,  p.  417,  1970 
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Since  the  TKE  monel  yields  the  turbulence  structure  (energy  and  length 
scale)  along  a streamline,  it  is  the  proper  model  for  the  study  of  free 
stream  turbulence.  In  the  Ttaci  and  Wilcox  paper,  the  stagnation  point 
problem  was  solved  by  dividing  the  flow  field  into  three  regions.  Region  1 
consists  of  a uniform  mean  flow  with  the  flow  unperturbed  by  the  presence 
of  a body.  The  turbulence  is  convocted  by  the  mean  flow  and  at  the  same 
time,  decays  by  dissipation.  Differential  equations  which  govern  this 
decay  process  are  obtained  from  the  conservation  equations  for  e and  X’. 

The  turbulence  level  obtained  near  the  body  is  used  as  the  free  stream 
turbulence  level.  In  Region  IT,  the  mean  flow  is  assumed  to  be  undisturbed 
by  the  turbulence  and  is  given  bv  the  lamina-  stagnation  point  inviscid 
solution.  I’he  turbulence  level  is  enhanced  by  the  mean  flow  strain  (even 
though  small)  in  this  region  as  shown  in  Figure  2.  (In  this  figure,  the  constant 
c IS  the  inviscid  flow  velocity  gradient).  The  solution  to  the  turbulence 
differential  equations  (note  the  mean  flow  is  assumed  to  be  known  in  these 
equa'-fons)  near  the  wall  will  be  the  turbulence  level  at  the  boundary  laver 
edge  which  is  the  input  valu«  to  Region  III.  Rogior  IH  is  the  viscous 
boundary  layer  region  where  the  flow  equations  are  completely  coupled  to  the 
turbulence  equations.  Solutions  for  this  region  yield  the  heat  tLansfer  to 
the  wall,  A numerical  method  of  solution  was  employed  in  lliis  region.  The 
match  or  patch  point  between  Region  II  and  III  was,  however,  somewhat  arbitrary. 

Solutions  for  both  the  cylinder  and  the  sphere  were  obtained  by  Traci  and 
Wilcox  as  shown  in  Figures  3 and  4,  Again  significant  heat  transfer  enhance- 
ment was  found,  the  amount  of  enhancement  depends  on  the  paramctc-rs 
where  D is  the  nose  diameter  and  Re^  is  the  Reynolds  number  based  on  I). 

In  sun.mary,  theoretical  studies  for  the  stagnation  point  heat  transfer 
enhancement  in  incompressible  flow  have  been  made.  The  theoretical  result.s 
agree  well  with  experimental  data  and  show  a significant  increase  in  heal 
transfer  due  to  free  stream  turbulence.  Corresponding  heat  transfer  theories 
for  S'.'pers''ni c flow  with  free  stream  turbulence  do  not  seem  to  be  available 
in  the  open  literature. 
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Section  3 


MASS  FLUX  AND  BULK  ENTHALPY 

The  heat  transfer  rates  measured  near  the  centerline  of  thi'.  RENT 
nozzle  are  roughly  2 to  2,5  times  higher  than  those  computed  using  the 
enthalpy  obtained  from  an  arc  chamber  heat  balance.  It  should  not  be 
Inferred  from  these  types  of  comparisons  that  the  stagnation  point  heat 
transfer  theory  predictions  are  low  by  a factor  of  2 to  2.5.  The  heat  rate 
distribution  across  the  nozzle  suggests  that  there  is  an  enthalpy  distri- 
bution across  the  nozzle.  In  order  to  determine  the  mass-averaged 
enthalpy  at  the  test  section,  the  enthalpy  flux  across  the  nozzle  must  be 
integrated,  A comparison  of  this  mass-averaged  enthalpy  with  the  enthalpy 
obtained  from  the  heat  balance  should  then  reveal  the  true  discrepancy 
betweex.  the  measured  and  predicted  heat  transfer  rate.  It  is  conceivable 
that  this  discrepancy  between  theory  and  data  is  not  a constant,  but 
varies  across  the  test  section.  It  does  not  appear  possible  to  determine 
this  fact  until  the  entire  problem  has  been  sorted  cut. 

The  mass  flux  and  the  mass-averaged  enthalpy  can  be  computed  if  the 
free  stream  total  pressure,  the  total  pressure  behind  the  shock  and  the 
corresponding  stagnation  point  heat  transfer  rates  are  known.  It  should 
be  noted  that  in  order  to  obtain  accurate  results,  the  profile  of  the 
total  pressure  behind  the  shock  should  match  the  heat  rate  profile.  Unfor- 
tunately, the  total  pressure  and  heat  rate  data  available  for  analysis 
during  this  program  were  not  perfectly  matched  and  the  results  presented 
below  must  be  considered  approximate. 

The  mass  flux  at  the  test  section  is  given  by 
r 

ra  = TT  ["  purdr 
e 

where  the  test  cross  section  is  considered  to  consist  of  two  semi-circles 
(-r^  to  0,  0 to  +1"^)  which  may  or  may  not  by  symmetric.  The  mass -averaged 
or  bulk  enthalpy  is  given  by 
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«B  = 


+r 

nf  ® pu  H rdr 
d -r 

e 


m 


The  above  expression  can  bo  rewritten  as 


_+r  , 

m = pM/0.604^  r dr 


+r  , , 

nj_^®  pMyo.604  ^ Hrdr 


«B  = 


m 


All  of  the  above  quantities  are  free  stream  values.  The  Mach  number, 
M,  is  obtained  from  the  ratio  of  total  pressures;  p,  the  static 
pressure  is  computed  as  a function  of  M and  the  free  stream  total 
pressure;  the  static  temperature,  T,  is  obtained  from  H,  the  local 
total  enthalpy,  and  Mi  the  compressibility,  Z,  is  obtained  from  the 
static  temperature  pressure.  The  local  total  enthalpy  is  determined 
from  the  Fay-Riddell  heat  transfer  relation 


q = CV_!2  (h  - h ) 

r w 

n 

It  is  this  last  relationship  which  requires  that  the  total  pressure 
(behind  the  shock)  profile  accurately  match  the  heat  rate  profile  in 
order  to  obtain  accurate  values  of  H. 

The  mass  flux  and  bulk  enthalpy  were  computed  using  the  above 
relations  and  the  total  pressure  and  heat  rate  data  repori.ed  by 
Brown-Edwards  (Ref.  3-1).  In  these  calculations  the  total  pressure 
behind  the  shock  was  assumed  j be  constant  across  the  test  section  at 
the  centerline  value  for  the  purpose  of  obtaining  the  total  enthalpy 
from  the  heat  rate  data.  This  approximation  was  necessary  since  the 
total  pressure  profiles  did  not  accurately  match  the  heat  rate  profiles 


Ref.  3-1  Biown-Edwards,  E.G.,  "Fluctuations  in  Heat  Flux  as  Observed 
in  the  Expanded  Flow  from  the  RENT  Facility  Arc  Heater", 
AFFDL  TR-73-102,  Nov.  1973 
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In  the  determination  of  the  Mach  number,  the  actual  total  pressure  pro- 
files were  used. 

The  mass  flux  and  bulk  enthalpy  computed  from  the  profile  data 
are  shown  in  Figure  5 for  three  nozzles.  The  approximation  discussed 
above  causes  the  mass  flux  to  be  high  and  the  bulk  enthalpy  to  be  low. 

It  is  difficult  to  determine  the  actual  error,  but  it  is  estimated  that 
the  computed  mass  flux  is  probably  about  10  to  20  percent  higher  than  the 
value  an  accurate  computation  would  yield  and  the  bulk  enthalpy  is 
probably  10  to  20  percent  too  low.  The  three  computation  points  shown 
for  the  1,11"  dia,  nozzle  correspond  to  the  average,  high  and  low 
level  heat  rates.  There  appears  to  be  no  significant  difference  between 
these  three  points.  .\lso  shown  in  Figure  5 are  the  bulx  enthalpies  derived 
from  arc  chamber  heat  balances  as  reported  b>  Brown-Edwards  (Ref,  3-1) 
and  Conley  (Ref,  3-2).  The  computed  bulk  enthalpies  a.e  about  10  and  30 
percent  higher  than  the  heat  balance  enthalpies  reported  by  Brown- 
Edwards  and  Conley,  respectively. 

It  appears  from  these  calculations  that  the  discrepancy  between 
the  measured  and  predicted  beat  rates  or  equivalently,  the  derived 
local  enthalpy,  may  be  somewhere  between  20  and  50  percent. 


Ref.  3-2  Conley,  H.  L, , "Enthalpy  Measurements  in  the  RENT  Facility 
Using  the  AEDC  Transient  Enthalpy  Probe",  AFFDL-TM-74-124- 
FXN,  May  1974 
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Section 

FREE  STREAM  TURBULENCE 


4,1  Overview 

This  section  considers  the  potential  effects  of  free  stream  turbulence 
on  stagnation  point  beat  transfer.  One  may  correctly  anticipate  that  such 
added  disturbances  will  increase  the  rate  of  heat  transfer  over  what  is 
predicted  by  a laminar  theory  like  Fay  & Riddell's  classic  formula  (Ref, 

4-1),  Although  there  have  been  no  direct  measurements  of  turbulence  in 
the  RENT  flow,  some  recent  work  by  Humphreys  (Ref.  4-2)  on  heat  transfer 
from  a rotated  arc,  which  also  has  added  swirl  because  of  the  mode  of  gas 
injection  into  the  arc  chamber,  shows  that  heat  transfer  to  a downstream  con- 
straining tubular  enclosure  is  increased.  Their  conclusion  was  that  "Large 
increases  in  the  heat  transfer  rates  to  the  tube  were  observed  with  rotation 
and  it  was  shown  that  these  increases  were  due  to  the  swirling  flow  and 

turbulence  caused  by  the  arc's  rotation".  The  Reynolds  number  of  the  RENT 

5 

jet  for  say  the  1,11"  nozzle  is  about  8 x 10  , an  extremely  high  value. 

Such  a jet  would  be  completely  turbulent  a short  distance  from  the  exit 
plane,  a process  which  could  be  further  shortened  by  the  known  arc  in- 
stabilities (Ref.  4-3)  or  the  turbulent  boundary  layer  that  probably  exists 
near  the  exit  plane  of  the  arc  chamber.  Because  the  model  diameter  may  be 
as  much  as  one-half  the  exit  nozzle  diameter,  fluctuations  in  the  jet  from 
any  of  the  flow  conditions  may  affect  the  boundary  layer  and  its  immediate 


Ref.  (4-1)  Fay  & Riddell,  "Stagnation  Point  Heat  Transfer  in  Dissc^iated 
Air",  J.A.o,  pp  73-85  (1958) 

Ref.  C4-2)  Humphreys,  J.  F.,  and  Lawton,  J,,  "Heat  Transfer  for  Plasma 
Systems  Using  Magnetically  Rotated  Arcs",  Journal  of  Heat 
Transfer,  p,  397-402,  Nov.  1972 

Ref,  (4-3)  Kesten,  J,,  and  Wood,  R. , "The  Influence  of  Turbulence  on 
Mass  Transfer  from  Cylinders",  Journal  of  Heat  Transfer, 
p.  321-327,  Nov.  1971 
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external  flow.  The  Reynolds  number  based  on  model  diameter  is  also  very 
high,  e.g,,  4 x 10^,  for  a quarter  inch  nose  radius  model  so  that,  as  will 
be  shown,  effects  of  external  stream  turbulence  would  be  accentuated.  It 
is  therefore  worthwhile  to  at  least  attempt  to  give  reasonable  bounds 
on  the  effect  of  turbulence  on  stagnation  point  heat  transfer  so  possible 
measurements  of  this  quantity  during  RENT  operations  can  be  quickly  evaluated. 

Our  discussion  will  first  center  on  the  extensive  experimental  results 
for  subsonic  flow,  then  go  on  to  the  limited  data  available  for  supersonic 
flow,  then  show  available  theoretical  results,  and  finally  make  recommenda- 
tions for  future  experimental  programs. 

4.2  Experimental  Results  for  Subsonic  Flow 


Experiments  on  effects  of  free  stream  turbulence  have  been  made  for 

cylinders,  spheres,  and  jets  impinging  against  plates.  One  of  the  most 

important  results  of  these  studies  is  that  the  effects  of  turbulence  are 

magnified  if  some  relevant  Reynolds  number  is  large.  If  the  turbulence 

level  is  expressed  by  the  definition  T = u'/U^,  where  u'  is  the 

fluctuating  level  of  free  stream  velocity  and  its  mean  value,  then 

increases  in  heat  transfer  relative  to  T = 0 are  expressed  in  terms  of 

1/2 

the  variable  T(Re)  or  T x Re  where  Re  is  the  Reynolds  number  based  on 
cylinder  or  sphere  diameter  or  discharge  nozzle  diameter.  Some  investi- 
gators have  used  mass  transfer  instead  of  heat  transfer  to  gage  the  effects 
of  turbulence  and  found  comparable  results. 

Kestin  & Wood  (Ref.  4-3)  have  summarized  results  for  a cylinder  as 

1/2 

shown  in  Figure  6,  The  value  of  Nu/R  at  T = 0 is  about  0.94  so  that 

e 

the  maximum  measured  effect  at  the  stagnation  point  is  (roughly)  a factor 
of  1,8  increase  in  the  Nusselt  number  hD/k.  A curve  fit  to  the  data  is 
shown  and  follows  (Eq.  4-1). 


.pj  . 0.945  + 3.48  \-S^/ 


TR  \ 

_S / 

1 nn  ' 


(4-1) 


This  expression  peaks  at  a value  of  TR^  of  43.6.  The  Reynolds  number 

for  a model  with  a half  inch  diameter  placed  in  the  exit  of  the  Mach  1.8 

5 1/2 

nozzle  has  an  R of  4 x 10  so  that  R = 600.  A 7%  free  stream  tur- 
e e 

bulence  level  could  then  produce  this  maximum  effect  on  the  stagnation  line 
of  a cylinder  in  cross  flow  if  the  effect  of  free  stream  turbulence  is 
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Figure  6 Measurements  for  the  Local  Rate  of  Heat  (Mass)  Transfer  at  the 
Stagnation  Line  of  Cylinders  in  Cross  Flow 
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independent  of  Mach  number. 
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An  important  feature  of  this  heat  (or  mass)  transfer  augmentation  is 
that  the  distributions  of  the  effect  around  the  stagnation  point  seem  to 
follow  a "displaced"  laminar  law.  This  is  illustrated  for  a cylinder  in 
Figure  7 from  Ref,  4-3,  Only  at  the  separation  point  at  6 = 105  degrees 
is  there  a difference  in  the  curves.  This  result  suggests  that  most  of 
the  boundary  layer  is  dominated  by  viscous  forces  and  that  turbulence 
represents  an  additional  "outer  flow"  transport  term  which  steepens  the 
relevant  wall  gradients. 

Early  work  by  Smith  & Kuethe  (Ref,  4-4)  shows  the  same  trends  and, 
for  clarity,  these  and  other  data  are  given  in  Figure  8,  One  of  the  key 
features  of  the  data  is  that  the  greatest  increment  in  heat  transfer 
experimentally  observed  was  a factor  of  1,7.  The  early  data  of  K.eston 
(indicated  by  the  O)  appears  to  jump  to  this  value  for  very  low  turbulence 
levels.  Therefore,  one  is  led  to  suspect  that  a maximum  value  may  indeed 
exist  and  that  the  simple  linear  trend  of  the  "theory"  is  not  correct. 

Results  for  the  augmentation  of  stagnation  point  heating  on  a sphere 
demonstrate  the  importance  of  having  the  turbulent  Reynolds  number  greater 

■3 

than  a value  of  about  7 x 10  . Since  the  RENT  model  Reynolds  numbers  are 
about  10^,  this  condition  could  be  met  if  free  stream  turbulence 
exceeded  TL, 

Some  results  from  Gostkowski  (Ref.  4-5)  are  shown  in  Figure  9 . The 
largest  measured  increment  of  heat  transfer  was  2, 5/1. 4 or  1.8  which  is 
in  apparent  agreement  with  the  results  found  for  cylinders.  The  straight 
line  showing  an  ever  increasing  effect  would  appear  to  be  a tenuous  extra- 
polation of  the  data. 


Ref.  (4-4)  Smith,  M.,  and  Kuethe,  A.,  "Effects  of  Turbulence  on  Laminar 
Skin  Friction  and  Heat  Transfer",  Physics  of  Fluids  No.  5, 
p.  2337-2344,  Dec.  1966 

Ref.  (4-5)  Gostkowski,  V.,  and  Costello,  F.,  "The  Effect  of  Free  Stream 
Turbulence  on  the  Heat  Transfer  from  the  Stagnation  Point  of 
a Sphere",  Int.  J.  Heat  and  Mass  Transfer,  ]_3>  PP-  1382- 
1386  (1970) 
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Figure  7 Effect  of  Turbulence  on  Mass-Transfer  Rate  A.round  a 
Cylinder  Reynolds  Number:  Re  = 125,000  (nominal) 
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Turbulent  Reynolds  number  Rej-Re  x Tg 


Fig«3?e 


Correlation  of  Free-Stream  Turbulence  Effects  on  Heat-Transfer 
to  Sphere 


Galloviay  and  Sage  (Ref.  4-6)  have  made  extensive  measurements  around  a 
spherical  body  as  well  as  at  the  stagnation  point.  Some  of  their  results 
are  given  in  Figures  10  and  11. 


The  Frossling  number  is  — . ; w so  that,  because  Nu  > > 2,  ratios 

7F  P ' 
e r 

between  the  quantity  yield  the  combined  effects  of  turbulence  and  Reynolds 
number.  The  quantity  a.  is  identical  to  T,  the  ratio  u'/U  . Notice  from 

L CO 

Figure  10  that,  at  the  largest  Reynolds  number,  68,000,  the  distribution 
of  Nusselt  number  for  all  turbulence  levels  follows  qualitatively  the 
distribution  for  the  1.3%  nominal  level  of  the  experimental  tunnel. 
Therefore,  the  increased  heat  transfer  level  at  the  stagnation  point 
followed  by  a monotomic  decline  is  valid  up  to  the  vicinity  of  separation. 
After  about  the  70°  angle  from  the  stagnation  point  there  are  differences 
in  behavior  depending  on  turbulence  level.  (It  should  be  noted  that 
these  and  other  investigators  have  cnecked  for  possible  effects  of  the  wake 
and  separated  flow  region  on  stagnation  point  heat  transfer  and  found  none. 
This  is  done,  for  example,  by  attaching  faired  streamline  shapes  to  the 
rear  of  the  body  which  suppress  the  recirculating  region).  The  maximum 
increase  in  stagnation  point  heat  transfer  was  1.55/1.15  = 1.35.  Looking 
at  Figure  il , the  maximum  change  in  heat  transfer  at  fixed  turbulence 
(cy^  - 0.256)  level  but  with  changing  Reynolds  number  was  found  to  be 
1.55/1.05  = 1.48.  There  seems  to  be  no  evidence  from  these  results  that  the 
effect  of  free  stream  turbulence  is  to  trip  the  laminar  boundary  layer  and 
make  it  turbulent.  Rather,  the  transport  in  the  usual  laminar  boundary 
layer  seems  to  be  augmented.  Aside  from  the  simple  "displaced"  form  of 
the  distribution,  this  conclusion  is  also  supported  by  the  Nu/,y^  scaling 
law  which  applies  even  with  the  added  free  stream  turbulence  level. 


These  increases  of  heat  transfer  appear  to  be  smaller  than  those 
encountered  for  the  cylindrical  case.  However,  this  is  not  necessarily 
accurate  as  shown  by  the  results  of  Ref.  (4-6)  given  in  Tabic  4-1. 


Ref.  (4-6)  Galloway,  T.R.,  and  Sage,  B.  II.,  "Local  and  Macroscopic  Thermal 
Transport  from  a Sphere  in  a Turbulent  Air  Stream",  A.I.Ch.E. 
Journal,  18^,  No.  2,  pp.  287-293  (March  1972) 
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Table  A-1:  Values  of  F = ^77! 

V^  P 
^ e r 


R 


T = 25.6% 


T = 1.3% 


68,000 

20,000 

5,237 


1.55 

1.38 

1.05 


1.19 

1.10 

0.9 


The  very  high  (1.37o)  level  of  tunnel  turbulence  has  itself  interacted  with 
the  variation  in  Rejmolds  number  to  produce  the  overall  variability  shown 
in  the  second  column.  If  the  low  value  of  0.9  for  F is  used  as  a reference 
then  the  increase  if  1.55/0.9  “ 1.72  which  for  the  = .256^68,000  = 

66.8  compares  with  cylinder  results  for  asymptotic  values  of  this  parameter. 
^■Jhen  the  low  value  of  0.9  is  used,  the  data  for  the  sphere  can  be  plotted 
as  shown  in  Figure  12.  The  low  point  has  a turbulent  Reynolds  number  of 
only  1340  which  apparently  confirms  the  results  given  in  Figure  9. 

The  situation  of  a free  jet  impinging  normal  to  a flat  surface  has 
been  studied  in  great  detail  by  Donaldson,  C.D.,  et  al..  Ref.  (4*7).  They 
show  the  ratio  of  heat  transfer  to  its  laminar  value  as  a function  of  dis- 
tance from  the  nozzle  exit  plane.  Since  the  jet  was  becoming  turbulent, 
the  results  represent  the  effects  of  turbulent  intensity  on  heat  transfer. 
They  are  shown  in  Figure  13,  where  R is  the  local  jet  Reynolds  number  = 
p^w^r^/p^.  The  subscript  c denotes  jet  centerline  values  and  r^  is  the 
jet  radius  at  the  impingement  location  corresponding  to  the  point  where 
the  jet  velocity  is  one-half  the  centerline  velocity.  Notice  that  the 
asymptotic  value,  corresponding  to  the  fully  developed  turbulent  jet  is  a 
ratio  of  2.2  which  is  slightly  higher  than  the  values  observed  for  spheres 
and  cylinders.  One  may  speculate  that  these  systematic  differences  are 
caused  by  the  more  severe  stagnation  point  pressure  gradient  encountered 
by  normal  flat  surfaces,  cylinders,  and  spheres  in  that  order.  When  these 


Ref.  (4-7)  Donaldsons,  C.D.,  et  al.,  "A  Study  of  Free  Jet  Impingement: 

Part  2 Free  Jet  Turbulent  Structure  and  Impingement  Heat  Trans- 
fer", Journal  of  Fluid  Mechanics  (1971)  Part  3,  pp.  477-512 
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SPHERES  (REF.  4-6) 
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flat  surface  results  are  cast  '.nto  the  form  used  for  cylinders  and  spheres 
the  results  are  as  shown  in  Figure  14.  The  general  trend  and  values  are 
in  agreement  with  those  found  previously  for  spheres  and  cylinders.  The 
scatter  of  data  indicates  that  the  parameter  is  insufficient  to 

correlate  results.  It  may  be  that  the  turbulence  scale  JL  ex;*ressed  in 
some  typical  nondimens ional  factor  like  (i/D)  may  also  be  important. 

If  the  RENT  facility  turbulence  level  is  connected  with  the  size  of  the 
exit  jet,  then  this  factor  would  be  important  when  comparing  heat  transfer 
results  at  the  different  Mach  numbers. 


A final  point  is  that  the  distribution  of  heat  transfer  coefficient 
follows  a laminar  type  behavior  away  from  the  stagnation  point  as  shown 
in  Figure  15.  Most  of  the  data  seem  to  be  shifted  by  a factor  of  about  2. 
The  authors  recommend  that  "for  most  engineering  computations  at  large 
Reynolds  numbers,  sufficient  accuracy  will  be  achieved  by  applying  the 
stagnation  point  turbulence  correction  factor  to  the  normal  laminar  heat 
transfer  rates  ...  it  will  be  found  that  the  heat  transfer  in  the  region 
near  the  stagnation  point  will  be  proportional  to  the  square  root  of 
Reynolds  number". 

Bearman  (Re'.  4-8)  shows  that  when  turbulence  scales  in  the  free  stream 
flow  are  very  mu;h  less  than  the  body  (model)  dimension  the  longitudinal 
component  of  turbulence  is  amplified  because  of  vortex  stretching.  This 
means  that  the  stagnation  region  of  a model  does  not  always  "damp  out" 
turbulent  external  flow  fluctuations. 


There  is  no  experimental  doubt  that  if  as  much  as  5%  turbulence  exists 
in  the  RENT  free  stream,  then  T^/^f«35,and  substantial  increases  in  stagna- 
tion point  heat  transfer  should  be  expected.  This  conclusion  is  not  changed 
oecause  of  comprtC<=ibility  or  the  presence  of  a bow  shock  wave  in  front  of 
the  model.  The  shock  wave  represents  an  adverse  pressure  gradient  to  the 
stagnation  region  streamline.  Some  directly  measured  turbulence  data  are 
given  in  Ref.  (4-9)  and  given  in  Figures  16  and  17. 


Ref.  (4-8)  Bearman,  P.  W.,  "Some  Measurements  of  the  Distortion  of  Tur- 
bulence Approaching  a Tv7c-Dimensional  Blunt  Body",  Journal  of 
Fluid  Mechanics,  pt.  3,  pp.  451-467  (1972) 

Ref.  (4-9)  Rose,  W.C.,  "Turbulence  Measurements  in  a Compressible 

Boundary  Layer,  AIAA  Journal,  12,  No.  8,  pp.  1060-1064  (1974) 
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For  these  particular  experimental  conditions,  the  pressure  rise  across  the 
shock  wave  was  2.5  times  the  free  stream  value.  The  static  pressure  rise 
in  the  RENT  facility  for  '■he  1.11",  Mach  1.8  nozzle  is  3.4  times  the  up- 
stream value  followed  by  a further  rise  of  a 1,25  factor  going  towards 
the  stagnation  point.  The  results  shown  in  Figure  17  indicate  that  the 
fluctuation  level  before  the  shock  wave/boundary  layer  interaction  was 
about  57.  and  after  the  shock  was  7 or  8%.  This  is  not  a spectacular 
change,  but  only  serves  to  shew  that  shock  waves  can  increase  at  least 
some  components  of  turbulent  flow.  Other  results  for  temperature  and 
density  fluctuations  may  not  be  directly  applicable  to  the  RENT  situation 
because  of  dissociation  and  other  real  gas  effects. 


Of  course,  detailed  turbulent  measurements  are  much  more  difficult  in 
very  high  enthalpy  streams,  but  the  effects  of  turbulence  of  heat  transfer 
can  be  inferred  from  some  reasonable  fluid  mechanical  aspects  of  these 
flows.  An  example  is  data  taken  from  the  inpingement  of  a "plasma"  torch 
onto  a large  vertical  plate  as  shown  in  Figure  18  from  Roma-Abu  (Ref. 
4-10). 


o PRESENT  STUDY  UPSTREAM  )C6.IOcm 
O PRESENT  STUDY  DOWNSTREAM  X=  9.65  cm 
0 GRANDE  upstream  X=  0.5  in 
C>  GRANDE  downstream  OF 
6”  SHOCK  X = 2.0  in 

MORKOVIN  DOWNSTREAM  OF 

EXPANSION  STATION  "F" 


MORKOVIN  UPSTREAM  STATION  ”A" 


— •■INCREASING  -i-  dp/dx 
O 


T7t  17  P" -P ♦*  P^oeciiT-rt 


Ref.  (4-10)  Roma'Abu,  M.M.,  "Heat  Pipe  Colorimetry  for  Plasma  Stagnation- 

Point  Heat  Transfer",  AIAA  Journal  10,  No.  3,  pp.  313-316  (1972) 
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PLASMA  BULK  TEMPERATURE 


These  heat  transfer  data  were  taken  with  a plasma  jet  of  argon  which 
operated  in  the  fully  laminar,  transitional,  and  fully  turbulent  jet  regime 
as  judged  from  the  Reynolds  number  at  the  jet  exit.  The  correlation  used 
was  that  developed  by  Incropera,  F. , and  Leppert,  G.,  (Ref.  4-11)  which 
related  transition  to  turbulence  with  large  increases  in  emitted  acoustic 
radiation.  The  top  two  figures  show  good  agreement  with  Fay  & Riddell 
theory.  As  the  flow  rate  is  increased,  transition  occurs  and  finally 
complete  turbulent  transition  in  the  jet  is  observed.  For  these  circum- 
stances, at  a plasma  bulk  temperature  of  5000°K  (approximately  one  of  the 
RENT  operating  values)  experimental  heat  transfer  values  are  1.6  times 
those  predicted  from  the  Fay  & Riddell  formula. 


Another  relevant  "hot  flow"  situation  is  that  described  by  O'Connor, 
et  al  (Ref.  4-12).  These  authors  use  a 1.5  megawatt  arc  heated  facility  to 
produce  a high  subsonic  flow  of  pure  nitrogen.  This  flui.d  is  allowed  to 
impinge  on  a vertical  plate  which  can  move  to  various  axial  locations.  At 
each  axial  location,  the  heat  transfer,  static  and  total  pressure,  and 
independently,  stream  enthalpy  were  measured.  These  quantities  were  then 
used  in  the  Fay  & Riddell  correlation  and  found  to  underpredict  the  observed 
heat  transfer  by  a factor  of  about  2.5.  These  results  are  shown  in  Figure  19 
where  the  stagnation  point  heat  transfer  is  plotted  versus  the  nondimens ional 
axial  distance  x/r.  The  solid  curve  uses  an  empirical  correction  factor 
(x/tj^)*^  in  the  Fay  & Riddell  correlation  type  of  equation,  viz: 


q . = 0.442 
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Ref.  (4-11)  Incropera,  F.,  and  Leppert  G.,  "Flow  Transition  Phenomena  for 

a Subsonic  Plasmajet",  AIAA  Journal  4,  No.  6,  pp.  1087-1088  (1966) 

Ref.  (4-12)  O'Connor,  T.  J.,  "Turbulent  Mixing  of  Axisymmetric  Jets  of 

r\  _•  ^ ^ ^ ^ ^,4  *.u  ^ 4 A 4 ^ 1 1 \\rnr\  ^ 

iaj.uxaxi.^  L/x oduCxci ecu  oxcxu^cii  wxuu  £i.itiL/xcue  rvxj.  ^ nvuu  icuii 

Report  RAD-TR-65-18  (1968)  AVCO  Corp.,  Wilmington,  Mass. 
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(r^  is  the  jet  radius  where  the  velocity  is  1/2  its  centerline  value.  Tt 
is  a function  of  the  axial  coordinate  x) . The  condition  of  total  enthalpy 
equal  to  4328  BTU  is  similar  to  RENT  operation.  The  Mach  number  0.782 
is  lower. 

4.3  Rough  Estimate  of  RENT  Turbulence  Level 

Because  of  the  demonstrated  importance  of  turbulence  (taken  together 
with  high  model  Reynolds  numbers)  it  is  appropriate  to  try  and  estimate 
the  turbulence  level  in  the  RENT  facility.  Only  very  limited  information 
on  fluctuations  is  available.  If  velocity  fluctuations  occur,  then  some 
fluctuations  in  total  pressure  would  result.  If  the  small  fluctuations 
given  by  Brown-Edwards  (Ref.  4-13  in  Figures  la,  lb,  Ic  are  viewed  to 
represent  such  fluctuations,  then  several  simple  models  may  be  used  to 
infer  the  shock  layer  turbulence  level.  A result  from  incompressible, 
homogeneous  turbulence  theory  (quoted  in  Ref.  4-8)  is  that 


.58 
pu 


If  p and  U are  evaluated  behind  the  shock  wave  and  ^ p^  = 13  psig  for 
the  1.11"  Mach  1.8  nozzle  (a  crude  estimate  by  inspection  of  the  figures 
in  Ref.  (4-13)  then  U'/U2  = .36,,  This  367o  level  which  appears  to  be  very 
high,  would  produce  the  maximum  (i.e.,  about  1.8)  effect  in  the  heat  trans- 
fer. Another  way  is  to  use  the  compressible  form  of  Bernouillis’  equation 
for  the  high  subsonic  flow  between  the  shock  and  stagnation  point.  A 
fluctuation  in  measured  total  pressure  would  then  relate  to  a velocit'y 
fluctuation  immediately  behind  the  shock  by  an  easily  derivable  relation 
ship; 


I 

i 


Ref.  (4-13) 


Brown-Edwards,  E.G.,  "Fluctuations  in  Keat  Flux  as  Observed 
in  the  Expanded  Flow  From  the  REN’i  Facility  Arc  Heater", 
Tech  Rep't  AFFDL  TR-73-102  (Nov.  1972)  A.F.  Flight  Dynamics 
Laboratory 
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Taking  the  static  pressure  behond  the  shock  to  be  6i.2  atm  (from  standard 
tables),  the  stagnation  pressure  to  be  76.5  atm  and  y “ 1.26  (its 

average  value)  we  obtain  for  the  1.11"  nozzle 

IT' 

=02  2. 

Uo  P 

2 '^o 

Since  Ap^  = (13  + 14.7)/14,7  = 1.88  atm  and  p^  = 76.5  atm  (slightly  lower 
than  actually  measured)  we  obtain  finally: 


= .054 


or  about  5.47,  turbulence  level,  a somewhat  more  realistic  figure.  This 
makes  the  paramet^’r  Ty^  = 37.8  which  would  give  a 507.  increase  in  heat 
transfer  rate. 

The  only  other  indications  of  fluctuating  flow  (aside  from  heat 
transfer)  are  given  by  Eader  (Ref.  4-14).  He  obtained  fluctuations  in 
the  intensity  of  copper  line  spectrometer  measurements.  No  power  spectral 
densities  of  the  signal  were  obtained  so  that  it  is  not  possible  to  say 
whether  the  flow  was  fully  turbulent,  i.e.,  a complete  distribution  of 
scales  of  motion  following  established  decay  laws. 

4.4  Theoretical  Predictions 

Theoretical  results  which  attempt  to  predict  the  large  observed  effects 
of  turbulence  on  laminar,  stagnation  point  heat  transfer  have  used  the 
simple  model  of  subsonic,  incompressible,  ideal  fluid  flow.  The  specific 


Ref.  (4-14)  Bader,  J.,  "Time  Resolved  Absolute  Intensity  Measurements  of 
the  5106  Angstrom  Copper  Atomic  Spectral  Line  in  the  AFFDLL 
RENT  Facility",  Tech  Report  AFFDL  TR-75-33  (Feb.  1975)  Air 
Force  Flight  Dynamics  Laboratory 
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form  of  the  turbulence  law  has  been  mixing  length  (Ref.  (4-4,  4-15)  or 
an  adaptation  of  turbulent  kinetic  energy  models  (Ref.  4-16).  In  all 
these  models  an  added  term  renresenting  the  eddy  diffusivity  of  heat  and 
momentum  is  used.  In  the  simple  ’’eddy"  law 

e = kU'y 
m 

where  k is  a constant,  U'  is  the  free  stream  fluctuation  level  and  y 

the  normal  coordinate.  In  the  more  complex  methods  e = e/uj  where  "e" 

m 

is  turbulent  energy  and  "uJ"  is  pseudovorticity  and  differential  equations 
must  be  solved  for  each  of  these  quantities.  The  more  complex  theory  has 
the  advantage  of  being  able  to  be  extended  to  compressible  flow  and,  more 
importantly,  to  account  for  the  correlation  lengths  in  the  flow.  Some 
typical  results  from  the  mixing  length  type  of  theory  are  shown  in 
Figure  20.  Figure  14  shows  the  theory  of  Ref.  (4-4). 

The  theory  continues  to  rise  with  Tv'^  whereas  the  experiments  reach 
an  asymptotic  value.  The  more  complex  theories  are  shown  in  Figure  21. 
They  appear  to  better  predict  the  trend  of  the  experimental  results.  Note 
that  the  correlation  length  appears  as  the  parameter 

Results  of  both  of  these  types  of  theories  for  spheres  give  about  the  same 
numerical  values,  although  the  experiments  tend  to  show  somewhat  lower 
influence  of  turbulence  on  heat  transfer  augmentation. 


Ref.  (4-15)  Galloway,  T.R.,  'Enhancement  of  Stagnation  Flow  Heat  and 

Mass  Transfer  Througu  Interactions  of  Free  Stream  Turbulence", 
AICHE  Journal,  19,  No.  3.,  pp.  608-617  (May  1973) 

Ref.  (4-16)  Traci,  R.M.,  and  Wilcox,  C.D.,  "Analytical  Study  of  Free- 

stream  Turbulence  Effects  on  Stagnation  Point  Flow  and  Heat 
Transfer",  AIAA  Paper  74-515,  7th  Fluid  & Plasma  Dynamics 
Conference,  Palo  Alto,  Calif,,  June  1974 
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Figure  20  Augmentation  of  Forward  Stagnation  Point  Transport, 

Comparison  with  Experiments  Involving  Circular  Cylinders  (Ref.  4-15) 
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Figure  21  Stagnation  Point  Heat  Transfer  to  a Cylinder  in  a 
Turbulent  Cross  Flow  (Ref.  4-16) 
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Another  approach  of  studying  the  effect  of  free  stream  turbulence  on 
the  stagnation  point  heat  transfer  is  the  so-called  vorticity  amplificat ion 
theory  advanced  by  Sutera  (Ref.  (4-17)  and  (4-18)).  This  theory  sugp-ests 
that  vorticity  amplification  due  to  stretching  of  vortex  filaments  in 
the  strongly  diverging  stagnation  point  flow  is  the  essential  mechanism 
underlying  the  strong  influence  of  free  stream  turbulence  on  heat 
transfer.  The  development  of  this  theory  was  based  on  the  study  of 
incompressible  two-dimensional  stagnation  point  flow.  In  this  study, 
the  free  stream  turbulence  was  mode''ed  by  a distributed,  unidirectional 
vorticity  in  the  form  of  a simple  sinusoidal  perturbation  of  a unique 
wave  number  k.  The  axis  of  this  added  vorticity  is  perpendicular  to  the 
surface  streamline  which  is,  according  to  the  vorticity  amplication  theory, 
the  only  direction  susceptible  to  ampli 'ication  by  stretching  in  the 
stagnation  point  flow. 

It  was  found  (Ref.  4-17)  that  amplification  depends  on  the  wavelength 
or  scale  of  the  added  vorticity.  Only  those  vorticity  components  with  a 
wave  length  \ greater  than  a certain  neutral  wavelength  (or  k < 1) 

amplify  and  can  therefore  affect  the  boundary  layer.  For  \ < or  k > 1, 
the  vorticity  is  dissipated,  for  > = or  k = 1,  the  vorticity  is  con- 
vected  toward  the  wall  with  no  net  amplification.  It  is  important  to 
emphasize  that  in  this  theory,  the  nature  of  the  vorticity  contained  in 
the  oncoming  flow  is  imposed  on  the  flow. 

Sutera  (Ref.  4-18)  has  made  calculations  for  the  stagna''i  m flow 
near  a cylinder  in  cross  flow  by  choosing  k = 1.5.  Increased  neat  trans- 
fer due  to  the  added  vorticity  was  found.  It  was  also  found  that  the  mean 
temperature  profile  is  much  more  responsive  to  the  added  vorticity  than  the 
mean  velocity  profile.  The  solutions  obtained,  however,  depend  on  a 


Ref.  (4-17)  S.  P.  Sutera,  T,  F.  Maeder,  and  J.  Kestin,  "On  the  Sensitivity 
of  Heat  Transfer  in  the  Stagnation  Point  Boundary  Layer  to 
Free-Stream  Vorticity",  J.  Fluid  Mech.,  Vol.  16,  part  4, 
p.  497,  1963 

Ref.  (4-18)  S.  P.  Sutera,  "Vorticity  Amplification  in  Stagnation  Point 
Flow  and  Its  Effect  on  Heat  Transfer",  J.  Fluid  Mech., 

Vol.  21,  part  3,  p.  513,  1965 
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so-called  amplification  parameter  A,  which  is  a free  parameter  not  provided 
b)’  the  theory. 

Based  on  Sutera's  work,  i^feeks  (Ref.  4-19)  has  applied  the  vorticity 
amplification  theory  to  study  the  effect  of  free  stream  turbulence  on 
hypersonic  stagnation  point  heat  transfer.  In  Weeks'  study,  the  dissipa- 
tion term  was  included  :n  the  <;nergy  equation,  the  wave  number  k was 
again  taken  to  be  1.5,  and  the  amplification  parameter  A was  related  to  the 
turbulent  intensity  by 

= IMhl)  M 

0.15  V 

o 

Where  is  the  mean  longitudinal  velocity  ahead  of  the  probe  shock  and 
|v|  is  the  velocity  fluctuation.  Based  on  the  Klebanoff  (Ref.  4-20)  values 
of  longitudinal  turbulence  intensity,  the  heat  transfer  through  a hyper- 
sonic nozzle  boundary  layer  was  calculated  and  compared  with  data  as 
shown  in  Figure  22.  It  is  seen  from  this  figure  that  the  calculated  value 
agrees  qualitatively  with  the  data  and  that  the  heat  transfer  with  free 
stream  turbulence  (q)  is  greater  than  that  without  free  stream  turbulence 
(q^).  Stagnation  point  heating  in  a plasma  jet  was  also  calculated  by 
Weeks  and  compared  with  data  as  shown  in  Figure  23.  Again,  qualitative 
agreement  was  found  between  the  measured  and  the  calculated  heat  transfer 
rates . 


Ref.  (4-19)  T.  M.  Weeks,  "Influence  of  Free  Stream  Turbulence  on  Hyper- 
sonic Stagnation  Zone  Heating",  Tech  Rept.  AFFDL-TR-67-195, 
May  1968 

Ref.  (4-20)  M.  V.  Lowson,  "Pressure  Fluctuations  in  Turbulent  Boundary 
Layers",  NASA  TN  D3156,  1965 


Figure  22  Variations  of  Stagnation  Point  Heat  Transfer  Through  a Turbulent 
Boundary  Layer 
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The  successful  application  of  vorticity  amplification  theory  to  a 
given  flow  problem  depends  largely  on  the  choice  of  wave  number  k 
and  the  relationship  of  the  parameter  A to  the  turbulent  intensity. 

The  theory  predicts  qualitatively  correct  results  as  shovm  by  Weeks. 

It  seems,  however,  difficult  to  obtain  quantitatively  accurate  results 
due  to  the  lack  of  any  universal  or  rigorous  method  in  determining  the 
two  parameters  k and  A.  It  should  be  pointed  out  that  the  vorticity 
amplification  mechanism,  which  is  certainly  one  of  the  important  mechanisms 
in  the  enhancement  of  heat  transfer  due  to  free  stream  turbulence  is 
included  in  some  of  the  modern  turbulence  models  such  as  the  turbulent 
kinetic  energy  model.  Such  a model  will  be  less  restrictive  in  its 
application  than  the  theory  of  vorticity  amplification. 


From  all  of  these  theoretical  and  experimental  results  it  is  clear 
that  any  free  stream  turbulence  in  the  exit  jot  of  the  RENT  facility  will 
have  (because  of  the  large  Reynolds  numbers)  a substantial  effect  on 
measured  values  of  stagnation  point  heat  transfer.  It  would  be  worth- 
while to  have  some  direct  measurement  (from  say  the  spectrum  of  fluctuating 
radiation)  of  its  existence  and  furthermore  to  extend  the  classical  heat 
transfer  results  of  Fay  & Riddell  to  include  effects  of  turbulence.  This 
should  be  done  using  a modern  formulation  of  turbulent-transport  which 
will  allow  inclusion  of  pertinent  (to  the  RENT  conditions)  real  gas  effects. 
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Section  5 

COPPER  PARTICLE  IMPINGEMENT 

In  order  to  determine  the  enhancement  of  stagnation  point  heat  trans- 
fer due  to  particle  impingement,  it  is  necessary  to  have  data  or  information 
on  particle  mass  loading  and  particle  size  distribution.  In  the  absence  of 
this  information,  only  an  upper  bound  estimate  can  be  made.  An  upper  bound 
estim.ate  was  made  based  on  the  following  two  basic  assumptions;  (a)  the 
injection  rate  of  copper  particles  from  the  electrodes  is  0.1  percent  by 
weight  of  the  gas  flow  rate  (Ref.  5-1);  (b)  the  number  density  of  gas 
phase  copper  in  the  model  shock  layer  is  10  part/cm  for  the  1.11  inch 

nozzle.  Ths  gas  density  in  the  model  shock  layer  for  the  1.11  inch  dia- 

-3  3 

meter  nozzle  is  approximately  2.8  x 10  gm/cm  , hence,  the  gas  phase 
copper  concentration  in  the  model  shock  layer  corresponds  to  0.017  per- 
cent by  weight.  For  other  size  nozzles  the  ratio  of  gas  phase  copper 
density  to  the  gas  density  in  the  shock  layer  is  assumed  to  be  the  same  as  the 
1.11  inch  nozzle. 

Let  us  denote  as  the  particle  density  in  front  of  the  shock, 

as  the  particle  impact  velocity  and  E^  as  the  collection  efficiency 
defined  as  the  percentage  of  that  can  reach  the  stagnation  region. 

For  unit  energy  accommodation  coefficient,  the  heat  transfer  from  particle 
impingement 


q.  = E D V/ 

imp  c c 1 /2 


(5-1) 


Both  E and  V.  are  functions  of  particle  size  (Ref.  5-2,  5-3).  Consider 

Cl 

first  the  case  of  uniform  particle  size.  By  assuming  an  approximate 


Ref.  5-1  "The  50  Megawatt  Facility"  by  Thermomechanic  Branch,  Flight 
Mechanics  Division,  Air  Force  Flight  Dynamics  Laboratory, 

Tech  Memorandum  TM-71-17  FXE,  Oct.  1971 
Ref.  5-2  G.  D.  Waldman  and  W.  G.  Reinecke,  "Particle  Trajectories, 

Heating,  and  Breakup  in  Hypersonic  Shock  Layers",  AIAA  J., 

Vol.  9,  No.  6,  June  1971 

Ref.  5-3  R.  F.  Probstein  and  F.  Fassio,  "Dusty  Hypersonic  Flows",  AIAA  J. 
Vol.  8,  No.  4,  April  1970 
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heating  and  evaporation  process,  together  with  the  known  injection  rate 
and  shock  layer  copper  number  density,  this  size  can  be  determined  by 
the  conservation  of  copper  mass  as  described  below. 

In  the  arc  heater,  we  assume  that  the  particles  are  spheres,  that 
the  particle  material  immediately  reaches  its  boiling  temperature  as  it 
comes  out  from  the  electrode,  absorbs  its  latent  heat  and  evaporates 
without  altering  the  local  heat  transfer  rate,  and  that  they  retain 
their  spherical  shape  as  they  evaporate.  With  these  assumptions,  net 
energy  conservation  for  a particle  subject  to  an  average  heat  transfer 
rate  q gives 

^dr 

q = - ( 

where  o is  the  particle  material  density,  r is  the  particle  radius, 
and  Q = latent  heat  of  sublimation  for  solid  particle  and  Q = latent 
heat  of  evaporation  for  liquid  particle. 


For  copper 


(oQ)  liquid 


1.12.  Hence  in  the  determination  of  par- 


ticle size,  there  is  little  difference  between  solid  and  liquid  particles. 
Solid  particles  in  the  heater  section  will  result  in  slightly  larger 
particles  in  the  shock  layer  and  therefore  will  have  a higher  impact 
velocity  and  higher  collection  efficiency.  For  the  purpose  of  estimating 
the  upper  bound,  we  will  consider  solid  particles  only. 

Assuming  that  the  particle  velocity  is  in  equilibrium  with  the  gas 
velocity,  or  equivalently,  that  the  Reynolds  number  based  on  the  particle- 
gas  slip  velocity  is  small,  then  from  Ref.  (5-4)  we  have  the  Nusselt 
number  - 2.  From  the  definition  of  Nusselt  number,  Eq.  (5-2)  can  be 
written  as 


Ref.  (5-4)  T.  Yuge,  "Experiments  on  Heat  Transfer  from  Spheres  Including 
Combined  Natural  and  Forced  Convection",  Journal  of  Heat 
Transfer,  Aug.  1960 


d(— ) 
'r 

o 


4k  (T  -T  ) 
o b 


(5-3) 


where  r^  is  the  radius  of  the  particle  ejected  from  the  electrode  (at 
t = 0)  j k is  the  heat  conduct' -'i  ty  of  air,  is  the  average  core  tem- 
perature in  the  arc  heater,  is  the  copper  boiling  temperature.  Assume 

- 6000°K,  which  is  approximately  the  temperature  in  the  gas  cap  over  the 
blunt  model.  At  6000°K,  k = 1.516  x 10  ^ . The  boiling  temperature 

for  copper  at  100  atm  pressure  is  4240°K.  ^¥*or  copper,  a = 8.9  gr/cm^  and 
Q = 8.14  X 10^  cal/mole.  Using  these  numbers,  Eq.  (5-3)  can  be  v.riiten  as 


(^)=  2.15  X 10-^ 

O r 


(5-4) 


The  residence  time  of  the  particle  in  the  heater  is  estimated  as  X/U^, 

where  U is  the  average  gas  velocity  in  the  heater  and  X is  taken  to  be 
o 

the  length  of  the  subsonic  portion  of  the  nozxle.  This  time  is  estimated 
to  be  one  msec.  From  Eq.  (5-4),  the  particle  size  in  front  of  the  model  shock, 


r , IS 
e’ 


r = r . 
e o 


yi  - 2.15  X 10"'^/r„2 


(r  in  cm) 
o 


(5-5) 


The  density  of  copper  vapor  immediately  behind  the  shock  due  to  the 
evaporation  of  copper  in  the  arc  heater  is 

-7  3/2 

^u(g)  = ^s  J 


(5-6) 


where  p is  the  gas  density  immediately  behind  the  shock. 

®o 

Let  US  consider,  for  the  time  being,  that  heating  in  the  shock 

layer  is  negligible.  In  this  case,  we  have  from  the  assumption  (b) 

-3 

mentioned  before  that  p . . - 0.17  x 10  p 

cu(g)  s^ 

Substituting  Eq.  (5-7)  into  Eq.  (5-6),  r^  can  be  determined.  Con- 
sequently from  Eq.  (5-5),  we  obtain  the  particle  size  in  front  of  the 


shock  as  r 


(5-7) 


Lz.o  pm. 


The  particle  density  before  the  shock  is  therefore 


p = 0.83  X 10 
c 


(5-8) 
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where  P^  is  the  gas  density  before  the  shock.  Equation  (5-8)  is  assumed 
to  be  applicable  to  all  nozzle  sizes. 

In  order  to  make  certain  that  heating  of  the  particle  in  the  shock 
layer  is  indeed  negligible,  a close  examination  into  the  heating 
mechanism  in  the  shock  layer  is  made.  The  heating  of  a particle  in  the 
shock  layer  is  governed  by  two  parameters  and  defined  as  (Ref.  5-2) 


H = EA/oQV  r 
e e e 

K = 3C  P A8or 
so  U so  e 


(5-9) 


E = 2(^) 
0 


M 


3.5 


[BTU/ft^^^-sec] 


where  V is  the  gas  velocity  and  M is  the  Mach  number  in  front  of  the 

® ® -3  3 

shock,  A is  the  shock  stand-off  distance,  = 2.498  x 10  slug/ ft 


is  the  drag  coefficient  of  particle  based  on  particle  cross- 


and 

sectional  area. 


Cj^  can  be  approximated  by  the  fo’.:mula  (Ref.  5-3) 


ae 


-n 


for 


Re  <'  1 


1 < R <10 
e 


a = 24 
3 


n = 1 
a = 24 


n ==  5 


(5-10) 


Re  > 10  a = 0.44  n 0 

where  Re  is  the  particle  Reynolds  number  based  on  the  particle-gas 
slip  velocity  and  the  particle  diameter. 

The  particle  impact  radius  is  plotted  as  a function  of 

and  K in  Figure  24  (from  Ref.  5-2).  The  collection  efficiency 
se 
m 

E (=  and  the  particle  impact  velocity  V.  are  also  plotted  in 
cm  1 

Figure!  25  and  26,  respectively  (from  Re^.  5-2).  In  Figure  25,  e is 
the  density  ratio  across  the  shock  e = — . 


,-3 


and  K is  of 
se 


For  a 12.8  lim  particle,  H is  of  the  order  of  10 
-2  ^ 

the  order  of  10  for  all  nozzle  conditions.  In  this  case,  we  see,  from 
Figures  24  and  26  that  not  only  is  the  heating  negligible,  but  also  that 
the  effects  of  the  shock  layer  on  particle  trajectory  are  negligible  small. 
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(i.e.,  E^  = 1,  V^  = V^).  In  fact,  shock  layer  heating  will  not  become 
important  until  the  particle  size  is  0.1  y,m  or  less.  Also  the  particle 
trajectory  will  not  be  altered  appreciably  in  the  shock  layer  until  the 
particle  size  is  1 ^m  or  less.  Based  on  these  considerations,  we  assume 


E = 1,  V.  = V (5-11) 

c re 

The  particle  density  and  the  particle  impact  velocity  are 

plotted  in  Figures  27  and  28,  respectively  for  various  nozzle  exit 

diameters  and  at  various  stream  enthalpys.  The  heat  transfer  due  to 

impingement  (from  Eqs.  (5-1),  (5-8),  and  (5-11))  is  presented  in 

Figure  29  along  with  the  range  of  obsezrved  convective  heat  transfer  4^. 

It  is  seen  that  within  the  estimated  range  of  the  stream  enthalpy  (from 

2000  to  6000^—),  the  impingement  heat  transfer  is  about  5%  of  the  ob- 
-lb 

served  convective  heat  transfer. 

The  results  of  Figure  29  represent  an  upper  bound  due  to  Eq.  (5-11). 
The  effect  of  the  particle  size  distrubution  and  the  uncertainty  in  the 
estimated  particle  residence  time  in  the  heater  can  only  reduce  both  the 
collection  efficiency  and  the  impact  velocity.  Therefore,  we  ccnclude 
that  the  effect  of  copper  particle  impingement  is  at  most  5 percent. 

This  conclusion  is  based  on  the  assumption  that  the  total  copper  mass 
fraction  is  0.1  percent  and  that  17  percent  or  0.017percent  by  mass  is  in 
the  vapor  state.  If  the  copper  mass  fraction  in  the  solid  p’.ase  is  higher, 
the  impingement  heat  transfer  would  be  proportionately  higher.  It  is 
therefore  desirable  to  obtain  data  on  copper  mass  loading  in  the  RENT 
flow  field. 


One  possible  method  for  obtaining  this  data  is  to  capture  the  copper 
particles  and  perhaps  the  copper  vapor  on  a membrane  filter  placed  in  an 
enthalpy  probe  at  a point  where  the  gas  temperatures  are  sufficiently 
low  to  enable  the  filter  to  survive.  The  copper  mass  captured  on  the 
filter  could  be  determined  by  emission  spectroscopy  or  atomic  absorption. 


of  spatial  enthalpy  gradients  on  stagnation  point  heat  transfer  were  per- 
formed to  assess  the  importance  of  these  effects  in  the  RENT  facility. 
Details  of  these  analyses  are  presented  in  Appendix  A and  B,  respectively. 
The  results  of  these  analyses  are  summarized  in  this  section. 

The  RENT  flow  field  contains  both  timewise  fluctuations  and  radial 


spatial  gradients  in  the  stream  enthalpy  as  demonstrated  by  the  heat  trans- 
fer measurements  results  shown  in  Figures  30  and  31  (Ref.  6-1).  The  heat 
transfer  results  shown  in  Figure  30  were  obtained  with  the  model  held  fixed 
near  the  nozzle  centerline.  Heat  transfer  models  are  normally  swept  across 
the  nozzle  to  prevent  destruction  by  excessive  heating.  The  results  shown 
in  Figure  30  were  obtained  when  the  traverse  carriage  malfunctioned  and  the 
model  stopped  for  about  40  msec  Heat  transfer  data  from  swept  models  nor- 
mally contain  the  effects  of  radial  gradients  as  well  as  timewise  fluctua- 
tions. When  the  radial  gradients  are  large,  as  is  the  case  in  Figure  31, 
the  timewise  fluctuations  can  be  seen  to  be  superimposed  on  the  radial 
gradient. 

In  Appendix  A,  the  effect  of  timewise  fluction  of  flow  properties  on 
stagnation-point  heat  transfer  is  investigated  by  considering  the  unsteady, 
one-dimensional  boundary  layer  energy  equation.  The  total  enthalpy  and 
free  stream  velocity  are  approximated  by  periodic  oscillations.  An 
analytic  solution  to  the  unsteady  energy  equation  is  obtained  by  the  method 
of  successive  approximations  up  to  first  order.  After  considerable  analysis 
and  manipulations,  the  final  result  for  the  increase  in  heat  transfer  due  to 
timewise  fluctions  is  developed  and  is  given  by  the  following  relation 


Ref.  (6-1)  Brown-Edwards , E.G.,  "Fluctuations  in  Heat  Flux  as  Observed 
in  the  Expanded  Flow  from  the  RENT  Facility  Arc  Heater", 
AFFDL-TR-73-102,  Nov.  1973 


m 


q„(t)  , 2 

= 1 - (1  - H^)  (q  ) r 

q (H  ) 
w e 


w 


Tn  the  above  expression; 


q^(t)  " time  averaged  heat  transfer  rate 

q^^(H^)  = heat  transfer  rate  corresponding  to  the  mean  total 
enthalpy 


q'  = heat  transfer  fluctuation  normalized  to  q (t) 
w w 

r = heating  enhancement  factor  (function  of  the  dimensionless 

frequency,  Bj.-) 

1 

In  the  Rent  facility,  the  heat  transfer  fluctuation,  q^  , is  of  the  order 
of  + 207o  at  a frequency  of  about  100  Hz  (Ref.  6-1).  For  these  conditions 
the  increase  in  heat  transfer  is  2%.  For  the  same  value  of  q^  and  a fre- 
quench  of  10^  Hz,  the  increase  in  heat  transfer  is  1.7%.  It  is  therefore 
concluded  that  enthalpy  fluctuations  in  the  RENT  facility  will  not  increase 
the  measured  stagnation  point  heat  transfer  results  by  more  than  a few  per- 
cent. 

An  analysis  of  the  effect  of  spatial  enthalpy  gradients  on  stagnation 
point  heat  transfer  was  also  carried  out.  This  analysis  is  presented  in 
Appendix  Bad  the  results  presented  below. 

By  considering  the  compressible  energy  equation,  it  is  shown  that  the 
increase  in  stagnation  point  heat  transfer  due  to  total  enthalpy  gradients 


Ws 


.' „ ,,  „ __  „ _ , „ .; 


where 


— = average  heat  flux  over  calorimeter  with  dH  /d0  ^ 0 

% 


— = average  heat  flux  over  calorimeter  with  dH  /d0  =0 
q s 

w 

0^  = hudy  angle  subtended  by  calorimeter  radius,  0(1°) 


dH 
s 

dC 


= total  enthalpy  slope  in  positive  0 direction 


= total  enthalpy  slope  in  negative 
o 


e 


direction 


The  above  equation  yields  some  interesting  results  on  the  effects  of  total 
enthalpy  gradients  on  stagnation  point  heat  transfer.  Note  first  of  all  that 
since  0 , the  body  angle  subtended  by  the  calorimeter,  is  generally  less 
than  0.1  radian,  the  change  in  enthalpy  will  have  to  be  rather  large  before 
there  is  a significant  increase  in  heat  flux.  (The  total  enthalpy  is  nor- 
malized to  the  value  corresponding  to  dH^/d0  = 0).  Secondly,  the  gradients 
must  be  of  opposite  sign  around  the  stagnation  point  or  the  effect  of  the 
gradients  will  cancel.  Thirdly,  it  is  possible  to  obtain  no  change,  an 
increase  or  a decrease  in  the  stagnation  point  heat  flux.  The  possibilities 
are  illustrated  in  Figure  32  , 

In  Figure  32a  the  total  enthalpy  gradient  is  a constant  so  that  the 
effect  of  the  gradients  cancel  according  to  the  above  equation.  In  this 
case  the  heat  flux,  when  total  enthalpy  gradients  are  present,  is  equal 
to  the  heat  flux  corresponding  to  the  average  total  enthalpy  across  the 
calorimeter.  For  the  case  shown  in  Figure  32b,  the  gradient  is  negative 
for  positive  0 and  vice  versa  for  negative  0 so  that  the  effect  is  an 
increase  in  heat  flux.  In  this  case,  the  actual  heat  flux  is  greater  than  the 
heat  flux  corresponding  to  the  total  enthalpy  averaged  over  the  calorimeter 
surface.  Figure  32c  shows  the  case  where  the  effect  of  total  enthalpy 
gradients  is  to  decrease  the  heac  flux. 
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BODY  ANGLE , 6 

(a) 


BODY  ANGLE,  9 
(b) 


BODY  ANGLE,  9 

(c) 


Figure  32  Possible  Effect  of  Total  Enthalpy  on  Stagnation 
Point  Heat  Flux 


^ - v'rt'M. 
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The  gradients  shown  in  Figure  31  can  increase  the  j/eak  heat  rate 
by  about  0.5%.  The  heat  rate  gradients  shown  in  Figure  31  are  typical 
of  RENT  data  and  it  is  therefore  concluded  that  spatial  enthalpy 
gradients  have  a negligible  effect  on  heat  transfer. 
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Section  7 
SWIRLING  FLOW 


In  the  present  context,  swirling  flow  refers  to  the  uniform  free  stream 
with  solid  body  rotation.  Nonuniform  flow  aspects  are  considered  in 
Section  4.  The  intent  here  is  to  investigate  the  effect  of  a swirl  com- 
ponent in  an  otherwise  uniform  flow.  An  analysis  of  the  effect  of  swirl 
has  been  performed  and  is  presented  in  Appendix  C.  The  results  of  this 
analysis  show  that  swirling  flow  effects  are  small.  Subsonic  heat  transfer 
data  to  a rotating  sphere  (Ref.  7-1)  were  located  which  confirm  the  results 
of  the  analysis.  These  data  are  shown  in  Figure  33  where 

RV  p 
R = — ^ 

^R 

RV  P 

R = 

e p 


V = R(ju 
r 

R = model  nose  radius 
(u  = flow  angular  velocity,  rad/sec 
= stream  velocity 

It  is  assumed  that  the  angular  velocity  in  the  free  stream  is  equal  to  the 

3 

angular  velocity  in  the  arc  chamber,  which  was  estimated  to  bo  10  rad/sec, 

_2 

the  ratio  of  R /R  is  of  the  order  of  10  . It  is  seen  from  Figure  31 

e e 

that  the  effect''^  of  circumferential  flow  on  heat  transfer  is  negligible. 


Ref.  7-1  Eastop,  T.D.,  "The  Influence  of  Rotation  on  the  Heat  Transfer 
from  a Sphere  to  an  Air  Stream",  Int.  J.  Heat  Mass  Transfer, 
Vol.  16,  pp.  1954-1957,  1973 
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Figure  33  Effect  of  Rotation  on  Heat  Transfer 


Section  8 

CHEMICAL  NONEQUILIBRIUM 

Based  on  our  review,  it  seems  certain  that  any  RENT  nonequilibrium 
effects  are  exceedingly  small.  There  would  at  most  by  only  a very  slight 
blocking  of  some  heat  transfer  resulting  in  an  underestimate  of  the  total 
enthalpy  from  Fay  & Riddell's  formula.  The  magnitude  of  the  effect  is 
less  than  5%. 

In  general,  several  requirements  must  be  met  for  a nonequilibrium 
blocking  effect  to  be  significant,  a)  A major  portion  of  the  stagnation 
point  enthalpy  must  be  associated  with  dissociation  or  ionization,  b) 
Homogeneous  chemical  times,  that  is,  times  required  for  the  atoms  or  ions 
to  recombine  in  the  gas  phase  must  be  greater  than  boundary  layer  flow 
times,  c)  Finally,  unless  there  is  considerable  ionization,  a noncataly- 
tic  surface  is  required  to  prevent  atom  recombination  on  the  wall  which 
otherwise  gives  a heating  rate  nearly  equivalent  to  the  equilibrium  case. 
With  ionization,  there  can  be  large  differences  in  heat  transfer  even 
with  a perfectly  catalytic  wall  (Ref.  8-1). 

Tables  8-1  and  8-2  provide  the  information  needed  to  address  Points  a) 
and  b) , Table  8-1  gives  a summary  of  equilibrium  stagnation  point  pro- 
perties for  "low",  "nominal",  and  "high"  values  of  RENT  total  enthalpy 
and  both  high  and  low  expansion  ratio  nozzles.  Flow  data  were  determined 
from  the  hypervelocity  air  flow  charts  of  Jorgensen  (Ref.  8-2),  compositions 


Ref.  8-1  Fay,  J.  A.,  and  N.  H.  Kemp,  "Theory  of  Stagnation-Point  Heat 
Transfer  in  a Partially  Ionized  Diatomic  Gas",  AIAA  J . , 

12,  2741-2751  (Dec.  1963) 


Ref.  8-2  Jorgensen,  L.  H.,  and  Baum,  G.  M.,  "Charts  for  Equilibrium  Flow 
Properties  of  Air  in  Hypervelocity  Nozzles",  NASA  TN  D-1333 
(Sept.  1962) 


are  from  the  Hilsenrath  (Ref.  8-3)  tables.  As  can  be  seen  from  Table  8-1, 
the  RENT  environment  is  fairly  mild  with  respect  to  a).  Ionization  is 
negligible  in  all  cases  and  dissociation  is  mostly  limited  to  oxygen  except 
in  the  high  enthalpy,  high  expansion  ratio  case  where  about  1/3  of  the 
nitrogen  is  in  atomic  form. 

Dissociative  kinetics  are  generally  fast  so  that  equilibrium  in  the 
shock  layer  is  a good  assumption;  however,  as  atoms  diffuse  into  the  cold 
BL  they  must  recombine  by  3-body  reactions.  If  the  gas  is  sufficiently 
rarefied,  the  recombination  process  will  be  slow,  thus  locking  up  some 
chemical  enthalpy  and  reducing  the  heat  transfer.  Table  8-2  lists  gas 
phase  Damkohler  numbers  (flow  time/chemical  time)  for  the  stagnation  con- 
ditions outside  the  BL  and  at  an  assumed  wall  temperature  of  1000°K.  Here 
the  flow  time  t^  is  the  time  for  a particle  in  the  free  stream  to  move 

one  nose  radius  r (Ref.  8-4). 
n 


t-  = r L 
f n 


2(P  - P )"^ 

S CO 


where  P and  p are  stagnation  pressure  and  density,  respectively,  and 
s s 

CO  refers  to  free  stream  conditions.  Since  oxygen  is  the  main  atomic 

species  the  chemical  time  t , is  the  characteristic  time  for  0-atoms 
‘ chem 

to  recombine  via 


0 + 0 -r  Mk  O2  + M 


and  is  given  by 


kX^CM]^ 

where  k = 3.8  x 10  ^*^T  ^exp(  ^^^)  cm^  molecules  ^sec  ^ 

= reaction  rate  constant 

r T . . .molecules, 

1_MJ  = total  concentration  ( x ) 


(Ref.  8-5) 


Ref.  (8-3)  Hilsenrath,  J.,  and  M.  Klein,  "Tables  of  Thermodynamic  Properties 
of  Air  in  Chemical  Equilibrium  Including  Second  Virial  Correc- 
tions from  i5Uo\  to  15,000"k,"  AEDC-TR-b5-58 , (March  1965) 

Ref.  (8-4)  Fay,  J.A.  and  Riddell,  F.R.,  "Theory  of  Stagnation  Point  Heat 

Transfer  in  Dissociated  Air",  J.  Aeronaut.  Sci.,  25,  121,  73-85, 
(1958) 

Ref.  (8-5)  Johnston,  H.S.,  "Gis  Phase  Reaction  Kinetics  of  Nuetral  Oxygen 
Species",  NSRDS-NBS20,  U.S.  Depart.  Commerce  (Sept.  1968) 
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X = mole  fraction  of  o-atoros  outside  B1 
o 


..3 


The  Datnkohler  number  F = tJt  , will  vary  as  1"'  in  the  BL  due  to 

f chem 

the  combined  effect  of  [m]  and  the  rate  constant  temperature  dependence. 

As  Table  8-2  indicates  with  high  expansion  ratio  nozzles  one  has  slow 

kinetics  (F  < 1)  in  the  free  stream,  but  still  fast  gas  kinetics  (F  > 1) 

at  the  cold  wall.  The  wall  region  is  most  important  in  determining 

heat  transfer  so  we  might  suspect  that  while  detailed  profiles  might  deviate 

from  equilibrium  in  the  outer  edgs  of  the  BL  the  heat  transfer  effects 

would  be  minimal  as  long  as  This  ir  borne  out  by  the  boundary 

layer  computations  of  Inger  (Ref.  8-6)  or  Fay  (Ref.  8-4)  which  show  almost 

-1 


no  difference  between  equilibrium  and  nonequilibrium  for  F^  > 10  *. 

Figure  34  taken  from  Inger  illustrates  this  and  shows  that  tc  achieve 
the  full  nonequilibrium  blocking  effect  required  F^  < 10 

Finally,  it  is  likely  that  even  the  small  (5%;  blocking  effects 

_2 

potentially  possible  according  to  Figure  34  for  F^  = 2 x 10  (high 
expansion,  high  enthalpy  case)  will  not  be  achieved  because  of  relatively 
high  surface  catalytic  efficiencies  (y^)  of  the  copper  models,  y^  has 
the  microscopic  interpretation  as  the  fraction  of  atoms  striking  che 
geometric  surface  of  the  model  that  undergo  heterogeneous  recombination. 
There  is  some  question  as  to  the  exact  value  of  y^  to  use  for  copper, 
but  available  data  indicate  it  is  in  the  range  of  0.1  " 1.0,  probably 
with  the  higher  values  applying  to  unoxidized  surfaces  (Wise,  Ref.  8-7; 


Ref.  (8-6)  Inger,  G.R.,  "Nonequilibrium  Stagnation  Point  Boundary  Layers 

with  Arbitrary  Surface  Catalycity",  AIAA  J.  Vol.  1,  No.  8, 
pp.  176  - 178,  Aug.  1973 

Ref.  (8-7)  Wise,  H.,  and  B.  J.  Wood,  "Reactive  Collisions  between  Gas 
and  Surface  Atoms",  Adv.  in  Atomic  and  Molecular  Physics. 

(D.  R.  Bates,  ed.,  291-353,  Academic  Press,  1967) 
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Figure  34  Nonequilibrium  Stagnation  Point  Heat  Transfer  Results  (Inger, 
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Carden,  Ref.  (8-8);  Anderson,  Ref.  (8-9);  Pope,  Ref.  (8-10). 

A Damkohler  number  for  the  heterogeneous  recombination  process  is 

expressed  as:  (Rosner,  Ref.  8-11) 


p \i  2(p  -p  j, 

2/3^  s s ■ s to).-|~-g 


where 


C = k p Sc 
w w w 


RT  1. 

I - / w .-S 

k - Y (^,r~) 

w w 2nM, 

A 

= surface  reaction  rate  constant 

7^  = surface  catalytic  efficiency  or  fraction  of  atoms  striking 

surface  which  will  recombine 

R = universal  gas  constant 

= wall  temperature 

M.  = mole  wt.  of  atoms 
A 

p = viscosity 
Sc  Schmidt  number  p/pD 

The  for  the  range  of  RENT  conditions  have  been  computed  and  are  also 

shown  in  Table  8-2;  as  can  be  seen  the  minimum  is  approximately  8. 

Refering  again  to  Figure  8-1  it  is  apparent  that  Q =8  would  behave 

w 

nearly  as  a perfectly  catalytic  wall  in  terms  of  heat  transfer. 

In  conclusion,  the  combination  of  relatively  modest  dissociation,  no 
ionization,  and  fast  gas  or  surface  reaction  kinetics  should  make  it  a cer- 
tainty that  there  will  be  no  significant  nonequilibri urn  effects  on  the  RENT 
heat  transfer  with  copper  models.  Even  with  models  v/ith  completely  non- 
catalytic  coatings  and  the  most  extre' e conditions  (H^  = 6000  B/lb, 
p^  =6.5  atm)  it  is  doubtful  if  a blocking  effect  could  be  detected. 


Ref.  (8-8)  Carden,  W.H.,  Heat  Transfer  in  Nonequilibrium  Dissociated 

Hypersonic  Flow  with  Surface  Catalysis  and  Second-Order  Effects", 
AIAA  J.,  4,  10,  1704-1711  (Oct.  1966) 

Ref.  (8-9)  Anderson,  L.A.  "Effect  of  Surface  Catalytic  Activity  on  Stagna- 
tion Heat-Transfer  Rates".  AIAA  J . . 11 . 5.  649-655  (May  1973) 

Ref.  (8-10)  Pope,  R.Pj.,  "Stagnation-Point  Convective  Heat  Transfer  ■’n 

Frozen  Boundary  T.ayers",  AIAA  J . , 6,  4,  619-626  (Apr.  1968) 

Ref.  (8-11)  Rosner,  D.E.,  "Scale  Effects  and  Correlations  in  Nonequilibrium 
Convective  Heat  Transfer",  AIAA  J . , 1,  /,  1550-1555  (July  1963) 
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Section  9 

TRANSPORT  PROPERTIES 


! 

i 


Tn  order  to  qualitatively  assess  the  importance  of  molecular  transport 
on  the  magnitude  of  stagnation  point  heat  transfer,  it  was  first  necessary 
to  establish  some  baseline  enthalpy  level.  Study  of  ReC.  9-1  shows  that  a 
level  of  4,500  BTU/lb  at  a chamber  pressure  of  100  atm  is  an  appropriate 
mean  upper  value.  For  these  ciroimstances  the  static  pressure  behind  the 
shock  wave  and  temperature  were  found  for  the  Mach  =1.8  nozzle  (1.11" 
diameter),  to  be  about  61  atm  and  524o\.  The  stagnation  conditions  are 
about  77  atm  and  5460°K.  For  these  conditions,  the  gas  consists  of  nitro- 
gen molecules  and  oxygen  atoms  and  some  oxygen  molecules  with  only  a small 
amount  of  atomic  nitrogen.  The  high  pressure  effectively  suppresses  the 
dissociation  of  nitrogen  molecules.  The  relevant  mole  fractions  of  the 
major  species  are: 


[01  = .261 
[N2I  = .692 
[02]  = .043 

For  these  equilibrium  stagnation  point  conditions  the  mole  fraction  of 
electrons  is  negligible  and  the  conductivity  of  the  gas  is  less  than  0.1 
mho/ cm. 

Heat  Transfer  Formulations: 

As  given  in  Ref.  9-2  the  standard  form  of  the  Fay  & Riddell  correlation 
for  heat  transfer  from  a simple  dissociated  gas  mixture  is  for  frozen  flow 


Ref.  9-1  Brown-Edwards,  E.G. , "Fluctuations  in  Heat  Flux  as  Observed  in 
the  Expanded  Flor  from  the  RENT  Facility  Arc  Heater',  AFFDL- 
TR-73-102,  Nov.  1973 

Ref.  9-2  Dorrance,  W.H.,  "Viscous  Hypersonic  Flow",  p.  302,  McGraw-Hill 
Book  Company,  New  York  (1962) 
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in  the  boundary  layer,  (catalytic  wall) 


« - I/O  j,.  1/2  (p  u )0.1 

(-q  ) = 0.763  Pr‘®-^(p  M [ (^)  ] 

x(H-h^)  [1  + i\P‘^^-\)  ^] 

The  only  difference  for  completely  equilibrium  flow  is  a very  small  change 
in  the  numerical  value  of  the  exponent  on  the  Lewis  number,  L,  from  0.63 
to  0.52. 

Effect  of  Non-Unity  Lewis  Number: 

Evidently  if  the  Lewis  number  is  unity,  the  added  contribution  to  the 
heat  transfer  from  the  degree  of  free  stream  dissociation,  expressed  as  the 
fraction  h^/H  of  total  enert,y  appearing  in  dissociated  gas,  is  eliminated. 
Based  on  the  calculated  shock  layer  composition  this  quantity  is  estimated 
to  be  0.241.  The  Lewis  number  is  a function  of  the  degree  of  dissociation 
as  shown  in  Figure  35  taken  from  Ref.  9-2.  The  value  L^  is  taken  to  be 
1.4.  For  the  given  stagnation  conditions  the  Lewis  number  is  then  found 
to  be  equal  to  1.25.  This  gives  only  a 3%  additional  heat  transfer  rate. 
The  largest  value  of  Lewis  numbar  is  1.464  corresponding  to  the  100  atm 
pressure  in  the  arc  chamber  with  dissociation.  This  value  gives  only 
a + 57o  increment  in  heat  transfer.  If  all  oxygen  were  dissociated,  the 
fraction  of  energy  in  this  form  would  rise  to  28%  of  the  total.  If  the 
largest  value  of  Lewis  number  is  used,  i.e,,  1.464,  a maximum  increase  of 
6.1%  would  occur.  These  conclusions  are  repeated  in  Table  9-1. 


Table  9-1  Effect  of  Non-Unity  Lewis  Number  of  Heat  Transfer 


Condi tion 


Equilibrium  at  Stag.  Point 
Best  Value  of  Lewis  Number 

Equilibrium  Stagnation  Point 
Largest  Value  of  Lewis  Number 

Frozen  Dissociation  Level  at  Reservoir 

Conditions  Largest  Value  of  Lewis  Number 


Heat  Transfer  Increase 


+6 . 1% 
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Efiect  of  Density-Viscosity  Product  Variations: 

A natural  quantity  that  occurs  during  analysis  of  boundary  layers  by 

means  of  similarity  variables  is  the  product  of  density  and  viscosity. 

Different  investigators  have  used  various  assumptions  about  how  this 

quantity  may  vary  across  a boundary  layer.  Sibulkin  (Ref.  9-3)  assumed 

pfi  = constant  (or  paT  so  that  the  ratio  is  unity  between  wall  and 

free  stream  conditions.  His  results  including  the  0.763  value  of  the 

constant  are  the  same  as  those  of  Fay  & Riddell  if  pu/pu  =1  and 

(D  O)  e e 

L =1.  Therefore  it  is  reasonable  to  assume  that  the  form  of  heat  transfer 
e 

correlation  is  not  affected  by  the  details  of  the  variability  of  pp  across 
the  boundary  layer. 

Fay  & Riddell  use  a Sutherland  type  viscosity-temperature  law  so  that 
po'T*^^  or  ppaT  ' \ The  term  ^ is  then  proportional  to 


T 

(f-) 

U) 


.024 


. The  largest  value  this  temperature  ratio  could  have  is  given  by 
the  ratio  of  shock  layer  temperature  to  ambient  temperature  which 
corresponds  to  the  heat  transfer  gage  not  rising  in  temperature.  This 
ratio  is  17:1  so  that  the  effect  of  neglect  of  this  term  is  to  underestimate 
heat  transfer  by  7%. 

Effects  of  Viscosity: 

The  "edge"  viscosity  appears  in  the  heat  transfer  equation  raised  to 
the  0.4  power.  Dorrance  (Ref.  9-2)  shows  typical  results  for  viscosity  as 
a function  of  temperature  as  reproduced  in  Figure  36.  It  can  be  argued  that 
one  should  use  the  value  of  viscosity  which  was  employed  to  derive  a particulai 
form  of  heat  transfer  correlation.  In  the  case  of  the  Fay  & Riddell  formula 
this  corresponds  to  Curve  (2)  of  Figure  36,  the  Sutherland  viscosity  law.  If, 
however,  (at  5000°K) , the  other  values  of  viscosity  were  to  be  used,  then 
relative  to  the  Fay  & Riddell  formulation  the  heating  rate  could  be  (pre- 
dicted) 167o  higher  or  9%  lower.  This  effect  is  essentially  the  variation 
of  "nose"  Reynolds  number  because  of  variations  in  viscosity.  Of  course, 
for  the  simple  equilibrium  conditions  taken  in  this  case,  the  density  is 
given  very  accurately  by  the  perfect  gas  law  applied  to  a dissociating  gas. 


Ref.  9-3  Sibulkin,  M.,  J.  Aeronautical  Sciences  1£,  No.  8,  pp.  570-571  (1952) 
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65 

Sutton  (Pef.  9-4)  treats  general  gas  mixtures  and  uses  a form  uryT* 
which  falls  within  the  stated  results. 

The  Prandtl  number,  which  appears  to  the  -0.6  power,  produces  no 
effects  worth  discussion,  since  the  thermal  conductivity  has  the  same  type 
of  temperature  dependence  as  viscosity.  It  varies  with  degree  of  dissocia- 
tion very  slightly  as  sho\m  in  Figure  35  with  a value  of  0.7  being  general- 
ly applicable. 

The  conclusion  can  be  drawn  that  for  the  assumed  conditions  of  equili- 
brium at  the  stagnation  point  the  variation  in  transport  properties  has 
very  little  effect  on  the  theoretical  heat  transfer  and  that  assuming 
Lewis  number  of  unity  has  no  serious  numerical  consequences.  These  conclu- 
sions would  have  to  be  altered  in  the  face  of  any  non-equilibrium  effects. 
For  example,  a large  excess  of  electrons  could,  if  their  temperature  were 
elevated,  produce  an  additional  heat  transfer. 


Ref.  9-^'  Sutton  K.,  and  Graves,  R.,  "A  General  Stagnation  Point 
Convective  Heating  Equation  for  Arbitrary  Gas  Mictures" 
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Section  10 
CONCLUSIONS 


The  results  of  the  studies  performed  during  this  program  lead  to  the 
following  conclusions: 

a)  The  discrepancy  between  the  bulk  enthalpy  obtained  from  an  arc 
chamber  heat  balance  and  that  computed  from  heat  transfer  and 
stagnation  pressure  measurements  appears  to  be  significant. 

b)  Free  stream  turbulence  can  account  for  this  discrepancy  if 
the  turbulence  level  in  the  RENT  facility  is  of  the  order 
of  a few  percent 

c)  Heat  transfer  enhancement  by  copper  particle  impingement  is 
negligible  if  the  copper  mass  loading  is  about  0.1%  and  of 
major  importance  for  mass  loadings  of  the  order  of  1% 

d)  Timewise  flow  property  fluctuations,  radial  enthalpy  gradients, 
swirling  flow,  chemical  nonequilibrium,  and  transport 
property  uncertainties  individually  have  a negligible  effect 

on  heat  transfer  in  the  RENT  facility.  Taken  in  total,  these 
phenomena  may  have  a noticeable  effect  on  heat  transfer. 
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APPENDIX  A 

EFFECT  OF  TIMEWISE  FLUCTUATIONS  OF  FLOW  VARIABLES  ON  STACNATION  HEAT  TRA.NSFER 
A.l  FORMULATION 

To  study  the  effect  of  timewise  fluctuations  of  flow  variables  on  stag- 
nation heat  transfer,  consider  the  unsteady,  one-dimensional  boundary-layer 
energy  equations,  which,  after  neglecting  the  viscous  heating  term,  can 
be  written  as 

M 4.  M iL_  /JJ  - i , 

bt  by  pby  Pr  by  p bt  (A--1) 


where  is  the  pressure  at  boundary  layer  edge,  H is  the  total 

enthalpy,  p is  the  density,  v is  the  velocity  in  the  y direction, 
which  is  the  outward  iiormal  to  the  surface,  Pr  is  the  Prandtl  number 
and  p,  is  the  viscosity.  In  order  to  simplify  Eq.  (A-1),  we  assume  Pr 
= constant,  p =■  p = p^,  and  that  v is  given  by  the  constant  den- 

sity inviscid  sol\iCion  for  sphere  near  the  wall,  i.e.. 


V 


wbeve  V^^  is;  the  free  stream  velocity,  R^,  is  the  nose  radius,  and 
f'  ^ 

— (density  ratic  across  a normal  shock} . 


In  this  iiwiscid  flow,  it  is  assumed  that  the  isentropic  conoxtion  is 
not  violated.  In  this  case,  we  have  from  the  second  law  of  thermodynamics 


<lhe  f 
etc 


d Pe 
, dt 


and  since  - < 1,  we  can  write  h " H 


+ 0 


(A- 2) 


Assuming  a perfect  gas,  we  have 


p h - H 

V e e y e e 


From  Eqs.  (A-3)  and  (A-4)  it  can  be  shown  that 
1 


p = C H ''-1 
^e  e ' 


where  C is  a constant,  and 


1 ^Pe_^«e 


p 5t  bt 


In  General,  and  can  be  represented  by  Fourier  Series  as 


V=(V)[1+S  ,Vi  cos  nft  + V'  sin  nft] 
~ ® o n = 1 In  2n 

CO 

CO 

H + (H  ) [1  + r , H'  cos  nft  + H'  sin  nft] 

e e o n = 1 


where  f is  the  normal  frequency. 


It  is  prohibitive  to  solve  Eq.  (A-1)  for  the  general  case  as 
described  by  Eq.  (A-7) , instead,  we  consider  the  simplified  case  of 
periodic  motion  as  given  by 


V =(V)  [l+V'  cos  ft] 

CO  00  o 


H = (H  ) [1  + H'  cos  ft] 

e e o 


Assuming 


, t,1/2  „ 1/2 

|i  ~ T ~ "e  > we  may  write 


H' 


p + (p,  ) [1  + ^ cos  ftl 

e e o ^2 


From  Eqs.  (A -5)  and  (A-8)  we  may  write 


H' 


p = (?  ) (1  + r cos  ft) 

'^e  e o Y-1 


Now  norr.ialize  H by  v by  (V^)^,  t by  and  y by 

Pr,  where  R is  the  Reynolds  number  defined  as 


(^)  (V  «N 

o o 


In  nondimensional  form,  Eq.  (A-1),  after  using  Eqs.  (A-6)  and 
(A-9)  becomes 


i'H  « ^rSe  /I  . I 


2 

= -S  ir  sin  t 


fp?  « 


where 


(A-10) 


(A-11) 


Assuming  H'  and  v'  to  be  small,  Sq.  (A-10)  can  be  solved  by  the  method 
of  successive  approximation.  To  illustrate  this  method,  we  rewrite 
Eq.  (A-10)  as 


Q M.  2 


|i,m.^,,^H.sint  = 2 ^yVcost-M 


, h:  (y-3)  ^ 

2 (Y-1)  ^^2 

The  zeroth  order  solution  is  found  by  neglecting  the  right  hand  side 
terms  of  Eq.  (A-12),  the  first  order  solution  is  obtained  by  evaluating 
the  right  hand  side  of  E^.  (A-12)  based  on  zeroth  order  solution. 
Similar  procedures  are  applied  for  higher  order  solutions.  We  will 
solve  Eq.  (A-12)  up  to  the  first  order.  Write 

H = H + H, 
o 1 

H satisfies  the  following  equation 
o 2 

OH  f-Tj — CiH  b H 

o ^ loe  o o _ „ I . , 

Yx  - ^ Vr  >'  X " ", 

t by 


(A-12) 


(A-13) 


(A-IA) 


A furtlipi-  approximation  is  made  to  simplify  the  equation  by  replacing 
the  coefficient  in  front  of  the  second  term  of  Eq.  (A- 14)  by  unity. 

This  approximation  implies  the  velocity  in  the  boundary  layer  is  assumed 
to  be  a constant  evaluated  at  a distance  y^  given  by  2 
This  approximation  will  not  affect  the  conclusion  of  this  study,  but 
will  enable  us  to  obtain  an  analytic  solution. 


Equation  (A- 14)  therefore  becomes 


't  5t 


sin  t 


(A- 15) 


y = 0.  ; y 


H =1  il  ' os  t 
o — 


Similarly  is  governed  by  the  following  equations 

bH,  6H,  bH  ut  / 6^ 

^ 1 1 1 _ I . o . H'  (y-3)  o 

ft  — — ^ = v'  cos  t •—  + T"')'  cos  t 

bt  by  by  2 (Y-1)  ^^o 


(A-16) 


y = 0 •=  0 


= 0 


The  heat  transfer  toward  the  surface  is  given  by 

c Ct)  = 

^w''  ^ by 


In  our  nondimensional  form,  it  can  be  shown  that 

1 


by/  y -- 


v/here 


IS  norma 


lized  by/p  \(V  ) (H  ) . The  purpose  of  the  present 
\ e)o  o e o - 


study  is  to  obtain  the  time  average  of  the  unsteady  heat  transfer  q (t)  = 

w 

_ 2tt 

compare  it  with  the  steady  heat  transfer  based  on 
the  average  free  stream  enthalpy  q^(H^).  The  ratio  of  to  q^(Hg) 

will  determine  the  effect  of  time  fluctuations  on  heat  transfer. 


j.;„»^s  *'  ‘ . 

^ '•  .>»>»WWN<!»1»W»«i>»~»»^WWMm  IB<I||  I !■■«»*  *'«^.  ■'^1— r _ ' i r*' T 


A.  2 SOLUTIONS 

liq.  (A-li>)  is  linear,  it  is  convenient  to  adopt  the  complex 
notation  to  write  liq.  (A-1^)  in  the  form 


5H  511  b H 

3 - -ir ^ = i 8 11'  ^ 

^t  bt  by  .,2  ■ t 


(A-17) 


y = 0 H = H , y ->  oo  H = 1 + ll'e 
o w o 

The  solution  to  Eq.  (A-1^)  will  be  the  real  part  of  the  solution  to 
Eq.(A-17).  Write 

H = (H  ),  + (H  )„  e^*“  ; where  (H  ),  and  (H  ).  are  functions  of  y (A-18) 

ooloi  o ot 


Substituting  Eq.  (A-18)  into  Eq.  (A-17)yields 

b^(H  ).  b(U  \ 

°1_  + ^ = 0 

2 5y 

by  ^ 

y = 0,(Hq)^  = y-co  01^)^  = ! 

“t«0>2  - 2^  ■ 

•'  oy 


(A-19) 


(A-20) 


y = 0,  (H  )„  = 0,  y - 00  (H  )„  = II' 


The  solutions  to  Eqs.  (A-1'9)  and  (A-20)  are 

<"oh  ■ - "»>  "w 

(11  ) -II'  [1  - o-<''2  + ‘h)’'] 
o I 


Hence 


= (1  - H ) (1  - e’^)  + H +11'  fcos(t)  - e ^ cos  (t-k,y)]  (A-21) 

w w ■ i 


re  li  2 

ki  = sin  , k^  = hil  + r^  cos  , r^  = [1  + (AS^)  'J'"’  ’ 

0 = tan 

o ^ 

By  using  Eq.  (A-21)  and  after  some  manipulations,  Eq.  (A-16)  becomes 


bH.  bH  0 IL  jl«  ( 3 

o i -J  = e "(1  - H ) (v'  - 

^t  bt  by  2 w 2(y-l) 


) ~'^2^ 
-)  cos  + H ' e 


• [^^(l  + cos  2t)  + sin  2tJ 


tt.  »35>1 
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with 


‘"l  " 2(^  0^1^-k/)  + K^V]  cos  k^y  + [v«k^ -k^k^£^V-3)  3 3^^  ^^y] 

0^2  = %{[  k^k^-  -JY-J)  -v*k^]  cos  k^y  + C2"'(^-l)^/’^2^^  ■*■ 


We 


again  write  in  complex  notation  as 


h " ^^^1^0  (np2®^ 


Since  we  are  only  interested  in  the  time  average  heat  transfer,  it  is 
’>nly  necessary  to  determine  The  other  two  terms  make  no 

contribution  to  the  average  heat  transfer. 


The  governing  equation  for  (H, ) is 


l"o 


-k2y 


by  ...2  01  (y) 


by 


y = 0,  (Hp^  = 0,  y-  » , = 0 

The  solution  is 

(Hi)q  = c(l-e'^)-H'J^e"^dyr^e^^2"^^^o'^(y)dy 


where 


c-HV0i(pp  + 92(e,)W)^|)^ 
= % 


_ k^2[ki^+  k2^+  2k2(k2-l)]-k2[(k2-l)(ki^-k2^)] 

5 2 2 2 

2{k^+  k^  )[(k2-l)  + k^  ] 


+ (k^^-  k2^)(k2-l)  -2k^\ 


2[(k2-l)^+  k^^] 


From  Eqs.  (A-21)  and  (A-24)  we  obtain  the  average  heat  transfer 
ratio  as 


H'v'e,  + 02 (^) 


(A-23) 


(A-24) 


(A-25) 


It  is  interesting  to  find  that  0^^  turns  out  to  be  a constant  independent 
of  frequency. 


Equation  (A-26)  states  that  the  timewise  fluctuation  will  increase 

the  time  averaged  heat  transfer  to  the  wall,  the  amount  of  increase  is 

2 

proportional  to  the  fluctuation  square  as  H'v'  and  (H’)  . It  is  very 
desirable  to  relate  these  fluctuations  to  the  heat  flux  fluctuations. 
This  will  be  discussed  in  the  following. 

For  high  Mach  numbers,  velocity  fluctuations  are  proportional  to 
the  square  root  of  total  pressure  fluctuations.  Assuming  v’  small, 
then  v'  can  be  approximated  as 


where  is  the  total  pressure  fluctuation.  For  isentropic  flow  of 
a perfect  gas,  the  total  enthalpy  fluctuation  is  directly  related  to 
P^*,  i.e.,  P^'  + H',  therefore,  we  have 

v'  = %H'  (A-27) 


In  order  to  relate  H'  to  the  heat  flux  fluctuation  q^'  consider 
the  zeroth  order  solution.  From  Eq.  (A-21)  we  obtain 

5H_ 


k2H’ 


1 — T-)  _ ,1  (1"H  ) [1  + „ cos  t + ] 

w by  y = 0 ' w'^  1-H 

w 


Write 


By  comparing  Eq.  (A-28)  and  (A-29)  we  find 

(1-H  ) q' 


H’  = 


w w 


By  using  Eqs.  (A-27)  and  (A-30),  Eq.  (A-26)  finally  becomes 


^ = 1+  (l-H^)(qV  r 


where  T is  the  heating  enhancement  parameter  defined  as 


(A-28) 

(A-29) 


(A-30) 


(A-31) 
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r = 


1 + 

Y-l  2 


('A-32') 


4k 


r is  plotted  in  Figure  A-1  as  a function  of  dimensionless  frequency 

Pf 


A.  3 CONCLUSIONS 

It  has  been  shown  that  when  the  free  stream  enthalpy  and  velocity 
fluctuations  can  be  approximated  by  periodic  oscillations,  the  time 
averaged  stagnation  heat  transfer  is  increased  over  that  obtained  with- 
out oscillation.  The  amount  of  increase  is  a weak  function  of 
frequency  and  is  proportional  to  the  square  of  the  heat  flux  fluctuation. 

In  the  case  that  the  free  stream  fluctuations  are  highly  irregular, 
these  fluctuations  should  then  be  considered  as  turbulence. 

According  to  Ref.  A-1,  q'^  found  in  RENT  was  20%  at  a frequency  of 

100  Hz.  Under  this  condition,  we  find  that  = 2 x 10  ^ < < 1.  and  the 

effect  of  this  fluctuation  on  5 is  6.7%,,  which  is  small.  Fre- 

5 ^ 

quencies  as  high  as  10  Hz  were  observed  in  Ref.  A-2,  and  for  this 

frequency,  = 0.2.  Again  assume  a 20%  heat  flux  fluctuation  and 

the  effect  of  this  fluctuation  on  q^  is  found  to  be  6.4%.  Therefore, 

for  20%  heat  flux  fluctuation,  the  effect  on  q is  small. 

w 


Ref.  A-1  G.  G.  Brown-Edwards,  ’’Fluctuations  in  Heat  Flux  as  Observed 
in  the  Expanded  Flow  from  the  RENT  Facility  Arc  Heater”, 
AFFDL-TR-23-102,  Nov.  1973 

t 

Ref.  A-2  Danie]  M.  Parobek,  ’’AFFDL  Developments  in  the  H^  and  q 

Measurements  in  High  Energy  Arc  Flow”,  presented  40th  Semiannual 
Meeting  of  the  Supersonic  Tunnel  Association,  24-26  Sept.  1973 
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APPENDIX  B 

EFFECT  OF  TOTAL  ENTHALPY  VARIATION  ON  STAGNATION  POINT  HEAT  TRANSFER 


In  this  appendix,  we  will  demonstrate  the  effect  of  spatial  variation 
of  total  enthalpy  over  the  surface  of  the  heat  transfer  gauge  on  the  stag- 
nation point  heat  transfer.  We  assume  that  the  velocity  vield  is  not 
affected  by  the  spatial  nonuniformity  of  enthalpy.  Therefore,  we  will 
consider  only  the  energy  equation) 


The  dimensionless  shock  layer  energy  equation  for  a compressible  gas 
with  constant  Prandtl  number  written  in  a set  of  similarity  variables 
(4,ti)  is:  (Ref.  B-1) 

2 


^(fM)  + 2R  p = U ^(1-p  );^(f^  ^)  + R P ’^s  „ 

bri  bn,  ^ ^ s ^ r bn  bn  e irr;-  H= 


U 


(B-1) 


where 


f = 


U ’ 
s 


R == 
e 


p (0)U  a 
00 


-.!f£ 


is  the  tangential 


U is  the  tangential  velocity  in  the  shock  layer, 

velocity  immediately  behind  the  shock  and  p (0)  and  (j  (0)  are  the  density 

s s 

and  viscosity  immediately  behind  the  shock  at  the  stagnation  line, 
respectively.  The  total  enthalpy  H is  normalized  by  that  behind  the 
shock  at  the  stagnation  line.  A sketch  of  the  geometry  is  shown  in  Fig. 
B-1.  Due  to  the  normalization,  the  shock  position  is  at  n 1 and  r)  = 0 


Ref.  B-1  Y.  S.  Chou,  "Locally  Nonsimilar  Solutions  for  Radiating  Shock 

Layer  about  Smooth  Axisymmetric  Bodies",  NASA  CR-1989,  March  1972 
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is  the  probe  surface.  The  heat  transfer  gauge  extends  from  0^  to 

-e„  (e„  <<  1). 

o o 


AXIS  OF 
SYMMETRY 


Figure  B-1  Flow  Geometry 


Equation  (B-1)  can  be  formally  integrated,  subject  to  the  boundary 

conditions  that  atri=l,H=H  atri  = 0,H  = H,as 

o’’  w 


H = G^(n)  + c^G^Cn)  + 


where 


■ J'o  ^ Jo  i>  + VrTT  f ■ 


f 5“  dT]"]dTl' 


= C ? if  -'o  f^ 


and  the  constant  is  given  by 

H -H  -G,(l) 
r = s w P '' 

1 G^(l) 


It  can  be  shown  that  the  convective  heating  is  related  to  constant  C, 


\-<-wr'>  <=i 

e s 


where  q is  normalized  by  p , U , H , U and  r are  dimensionless 

W 00  00  00  s S 

quantities,  normalized  by  U , ...  i „ , . 

oD  and  a respectively..  Hence,  we  obtain 
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(B-5) 


dH 
s 

de 


0 C'^  _ 


dH 

g 


CjH  11-H^-G^(1)]G'2(1)H 


= 0 


dH 


Where  subscript  ( )’  denotes  the  solution  for  the  case  0. 


dH 


Note  for  the  constant  enthalpy  case  (i.e.,  — jr  - 0),  H - 1 due  to  the 

Qo  S 

normalizatign.  ii  is  the  average  H over  the  gauge. 

- _ -^o  ° H dm 
H *“ 

2^0 

Assuming  the  effect  of  free  stream  enthalpy  nonuniformity  on  the  velocity 
field  as  well  as  H^  can  be  neglected,  we  then  have 


G'^(l)  - G^(1)  , 0*2(1)  - G2(1) 


(B-6) 


For  small  0 we  may  write 
dP 

He  = 1 ® + ■ — 

s dO  0 


After  using  Eqs.  (B-6)  and  (B-7),  Eq.  (B-5)  becomes 


(B-7) 


(q„)  hh^ 


dH 


= [1  + 


s)0 

d0 


w 


d0 


=•  0 


= lo  f '-^Vr-f-  “df  lo  f d: 


r <%> 

s dH 

— - 0 
d0 

dr 


62(1) 


-]  T 
H 


']d'n 


(B-8) 


now 


Assuming  £ - 1 then,  ^2^^^  “ 2 


_ „ er£  ^ ^) 

"g 

I r dr 

Vr  p _s.  _s 

e r U d§ 

s ^ 

dr  dr 

for  R^  > > 1,  r^  ~ 0 , Ug  ~ sin  6 ~ e,  -j|  ~ ~ 1 

'Al  ^ 


Hence  for  R > > 1,  G,(l)  ~ 2 


^/k  P 
^ e r 


^ = 0 


- 0.65  + 0,330  where  (-)  sign 


From  Ref.  B-2,  we  obtain  V^g(<l^) 

is  for  9 > 0 and  (+)  sign  d0 

for  9 < 0.  Equation  (B-8)  is  now  reduced  to  (P^  = 0.72) 


dH 
s 

d0 


0 


/ N dH 


dH 

= ri  + d£-2 — T i 

0.65  + 0.330"  - 

" n 


(B-9) 


The  (••)  sign  in  the  denominator  of  Eq,  (B-9)  is  for  9 > 0,  (+)  sign 
is  for  0<O.  Let  us  write 

dH 


■•w 


dH 
s 

d0  ?t  0 


d0 


p®°  . N - nd0 

J-0o%^  d0  " ° 


200  ' *w  200 

then  q^  is  the  average  heat  flux  over  the  heat  transfer  gauge  in  the 


Ref.  B-2  R.  T.  Davis  and  I.  Flugge-Lotz,  "Second  Order  Boundary  Layer 
Effects  in  Hypersonic  Flow  Past  Axisymmetric  Blunt  Bodies" 

J.  Fluid  Mech.,  Vol.  10,  Pt.  4,  Dec.  1964 
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I 


Cv- 

fv: 


f 5;' 


clH 

case  that  — ^ 0 while  q^  is  the  average  heat  flux  over  the  gauge 


dll 

~ = 0 
d0 

From  Eq.  (B-9)  we  have 
0.48 


•^w  200  <1  -0o 


00 


d0  ^o  d0 

0.65  - 0.330 


- r,  , 0.24  p0o 


0o 


0 d0 

dO-’o - 

0.65  - 0.330- 

•T 


(B-10) 


After  0 integration,  we  finally  obtain  from  Eq.  (B-10) 


dH 


^ 1 + 0.18  0o[l  + O.520o]f-^)^] 

_ _ 


(B-11) 


‘w 

dH 


1 + 0.25  60  [-3f)^  ] 


where  denotes  the  slope  of  in  the  positive  0 direction, 

dli  - dHg 

— jTr)  indicates  — rr)  in  the  negative  direction.  Based  on  Eq.  (B-11) 
d0  o d8  o 

the  effect  of  the  total  enthalpy  spatial  variation  can  be  discussed  as 
the  following. 

a.  If  the  enthalpy  varies  linearly  over  the  heat  transfer  gauge  as 
shown  in  Figure  B-2a,  there  is  no  net  effect  on  the  heat  transfer. 

b.  If  the  gauge  is  located  on  an  enthalpy  maximum  as  shown  in 
Figure  B-2b,  the  effect  is  that  the  average  heat  transfer  in  this  case 
will  be  greater  than  that  of  the  constant  enthalpy  case. 


i 4 

t ,, 


The  amount  of  decrease  or  increase  from  the  constant  enthalpy  heat 
transfer  in  Cases  (b)  and  (c)  obviously  depends  on  the  local  slope 
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APPENDIX  C 

EFFECT  OF  SWIRLING  FLOW  ON  STAGNATION  POINT  HEAT  TRANSFER 


The  effect  of  free  stream  swirling  motion  on  stagnation  point  heat 
transfer  to  a blunt  body  is  considered  in  this  Appendix.  It  is  assumed 
that  the  effect  of  coupling  between  compressibility  and  swirling  motion 
on  heat  transfer  is  negligibly  small  so  that  an  incompressible  analysis 
will  give  the  correct  relative  change  in  heat  transfer  due  to  swirl.  This 
analysis  is  therefore  for  a viscous  incompressible  fluid  with  a swirling 
motion  superimposed  on  a uniform  flow. 

For  the  incompressible  swirling  motion  near  an  axisymmetric  stagna- 
tion point,  the  inviscid  velocity  components  can  be  written  as 


U = Ar  , V = Or,  W = 2AZ 

where  U is  the  component  in  the  r direction,  Figure  C-1,  v in 
the  cp  direction,  W is  the  component  in  the  negative  Z direction, 


(C-11 


Figure  C-1  Coordinate  System 


n is  the  angular  velocity  of  a solid  body,  and  A is  a parameter 
related  to  the  strearawise  velocity  gradient  at  the  stagnation  point. 


A = j.  - Q*  plane  Z = 0 is  the  stagnation  plane. 

Z = 0 

From  Eq.  (C-1)  and  the  inviscid  momentum  equation,  we  obtain  the 
pressure  gradient  as 

■ p S " (a"-  nSr  , - ^ 7^  = 2Afir  , - 0 (Boundary  Layer) 

^ (approximation) 

The  boundary  layer  equations  in  this  case  can  be  written  as 
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r 


N.«  V. 


(A^  . fiS  r + V [^  + 1^  6 + ^] 
5r  bZ“ 

U + w ^ + 2Aflr  + u r-^  + 7-  (-)  + ^ ! 

'br"  ^ bZ- 


uM+w^-^ 

br  bZ  r 


5r  " bZ 


f + - +r^  = o 

br  r bZ 


The  boundary  conditions  are 


(C-2) 

^C-J) 

(C-^t  ' 


Z-0  , U-0  ,W-0  (no  slip  condition) 

Z-00  , U-Ar  , v=Or 

Define  the  dimensionless  variable  Q ~ 


(C-5) 


and  assume  the  velocity  components  in  the  boundary  layer  having  the  form 
U = r (A  + Q)  F(0  , V = rnC(0  , W = -27(A+r)u  H(Q  (C-6) 

We  obtain  equations  for  H,  G,  and  F from  Eqs.  (C-2)  - (C-4)  as 


H'"  + 2HH"  - CH')^  + = 0 

(A«)2 


2A 

G"  + 2HG'  - 2H'G  + ^ = 0 


F = H' 

Where  the  subscript  ( )'  denotes  differentiation  with  respect  to  Q. 
The  boundc 

; = c 

c 

The  solution  to  Eqs.  (C-T)  - (C-10)  yields  the  velocity  field  for  a 
swirling  stagnation  point  boundary  layer. 

The  stagnation  point  energy  equation  can  be  written  as 

g"  + 2P^Hg'  = 0 
T - T 

w 

where  g = — temperature  at  wall  and  '\\  the 

"st  ' ^W  ” 

stagnation  temperature  behind  the  normal  shock. 


conditions 

for 

H and  G 

are 

H = 0 , 

H' 

= 0 , G 

^0 

H = 

AMI 

G 

= 1 

(C-?) 

(C-6) 

(C-9) 


(C-10; 
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Til--  I oiinclary  conditions  for  g are:  Q = 0 g = 0 

1.  Now  the  heat  transfer  to  the  wall  is 

V “'S  = ’"’’st 

V 


Hence 


-^-T — b'(0) 

n-o  ~ 


(C-12) 


Since  the  energy  equation,  Eq.  (C-11),  is  coupled  with  Eqs,  (C-7)  - 

(C-10)  through  the  velocity  function  H,  g'(0)  will  be,  in  general, 

a function  of  ~ 

Q ^ 

For  small  — we  see  that  Eq.  (C-7)  is  decoupled  from  Eq.  (C-8)  to  the 

^ n 2 

order  of  (— ) . Also  the  boundary  condition  for  H as  ^ » can  be 

approximated  as  H'  - 1.  in  this  case,  the  equation  for  H becomes 


H' 


+ 2H  H' 


(H')  +1=0 


(C-ll- 


# 

T. 

X, 

tr 


C = 0 , H = 0 , H'  = 1 

Equation  (C-13)  is  identical  to  the  equation  for  the  non-swirling  case. 
The  solution  to  Eq.  (C-13)  is  available  in  Ref.  (C-1),  Hence,  for  small 

we  have  H - + O(^)  which  implies  g'(0)  - g'(0)^_Q+  O(^)  and 

the  effect  of  sv^irl  on  heat  transfer  for  small  — is 

A 


= yHO[l  + 0(J)J  J<<1  (C-14) 

^'lw^n=0  A A A 

From  Ref.  (C-2)  the  parameter  A is  found  to  be  related  to  the  density 
across  tlie  normal  shock  as 


A = (3e)  05  where  e 


We  finally  obtain 


R - nose  radius 
s 


% 


[1  + 0^)]  for  I < < 1 


(C-15) 
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Jvf  ^<,?..sy  ;w  ’ 


From  this  analysis  wo  soo  that,  physically,  the  effect  of  swirling  motion 
on  stagnation  poont  heat  transfer  is  two-fold.  The  first  is  to  increase 
the  boundary  layer  thickness  duo  to  the  influence  of  the  (positive 
direction)  swirling  secondary  flow.  Tliis  in  turn  will  decrease  the 
heat  transfer.  The  second  effect  is  to  increase  the  streaira^/ise  pressure 
gradient  at  the  boundary  layer  edge  due  to  the  centrifugal  force.  This 
v;ill  increase  the  heat  transfer.  The  second  effect  dominates  so  that  the 
net  result  of  these  two  efiects  is  to  increase  the  heat  transfer  at  the 
stagnation  point.  This  increase  in  heat  transfer  is  small  for  the  RENT 
flow  conditions  as  discusser  previously. 


Ref.  C-1  H.  Schlichting,  "Boundary  Layer  Theory",  4th  edition,  McGraw- 
Hill,  1960,  pp.  81-83  and  pp.  177-179 

Ref.  C-2  Y.  S.  Chou,  "Inviscid  Hypersonic  Flow  Past  Blunt  Bodies",  AIAA  J. 
Vol.  7,  No.  1,  Jan.  1967 
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